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Abstract:  Acoustic  propagation  and  diffraction  of  high-amplitude,  short 
duration  signals  through  a  forest  edge  have  implications  for  noise 
mitigation  and  battlefield  acoustic  sensors.  While  the  acoustic  significance 
of  this  unique  environment  has  been  noted  in  the  past,  it  has  not  been 
studied  in  any  detail.  Acoustic  signals  that  have  propagated  through  a 
forest  edge  yield  complicated  pressure-time  histories  for  receivers  both 
within  and  outside  the  forest.  Several  physical  processes  contribute  to  this 
complexity,  including  the  physical  structures  of  the  biomass  and  ground 
and  the  microclimate.  A  deep  understanding  of  acoustic  propagation 
through  this  unique  environment  may  lead  to  strategic  placement  of  fire 
breaks  for  noise  mitigation  and  improved  signal  processing  algorithms  for 
use  with  acoustic  detection,  and  direction-finding  and  range-finding 
sensors.  Because  of  the  broad  scope  of  issues  that  could  be  addressed  once 
acoustic  propagation  and  diffraction  at  a  forest  edge  is  understood,  it  is 
important  to  study  this  unique  environment  in  detail.  This  report  provides 
documentation  of  a  field  experiment  conducted  as  part  of  a  study  of  the 
acoustic  properties  of  the  forest  edge  environment. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 

DESTROY  THIS  REPORT  WHEN  NO  LONGER  NEEDED.  DO  NOT  RETURN  IT  TO  THE  ORIGINATOR. 
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1  Introduction 

The  experiment  described  in  this  data  report  took  place  on  Range  79  at 
Camp  Ripley,  MN  during  19-27  June  2006.  Data  recordings  were  ob¬ 
tained  on  22-25  June  2006.  The  location  was  chosen  because  of  its  many 
favorable  attributes.  It  is  a  uniform  pine  stand  with  a  distinct,  straight 
edge.  The  ground  is  flat  with  a  sufficient  open  field  adjacent  to  the  stand. 
The  site  was  easily  accessible  and  of  sufficient  size. 

Objectives 

This  experiment  had  several  goals,  as  stated  below: 

1.  Obtain  acoustic  measurements  in  the  open  field,  within  the  forest  edge 
area,  and  within  deep  forest.  Measure  at  more  than  one  source  and  re¬ 
ceiver  height  combination.  Use  these  measurements  to  track  the  acoustic 
propagation  through  these  three  distinct  media. 

2.  Determine  the  changes  in  ground  impedance  characteristics  and  the  effect 
on  acoustic  propagation  as  sound  travels  through  the  forest  edge. 

3.  Obtain  measurements  of  the  acoustic  backscatter  from  the  forest  edge  into 
the  open  field. 

4.  Investigate  the  changes  in  meteorological  profiles  as  one  moves  from  the 
open  field,  through  the  forest  edge,  and  into  the  forest.  Take  sufficient 
measurements  to  verily  the  theory  and  computer  model  calculations  of 
Arnold  Tunick  (Army  Research  Laboratory  meteorologist). 

5.  Determine  the  importance  of  the  biomass  structure  at  the  forest  edge  on 
acoustic  propagation. 

6.  Determine  a  characteristic  size  for  the  effective  transition  zone  (i.e.,  the 
distance  needed  to  transition  from  open  field  to  forest  conditions  and  vice 
versa). 

Approach 

Four  acoustic  sources  were  used  during  the  course  of  this  experiment: 

1.  1.25-lb  bricks  of  Composition  C-4, 

2.  a  propane  cannon, 

3.  a  loudspeaker  using  a  white  noise  signal,  and 

4.  a  .45  cal  pistol  with  blanks. 
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These  sources  were  chosen  to  represent  a  broad  variety  of  source  types 
(impulsive  and  continuous,  linear  and  nonlinear)  and  the  full  range  of  fre¬ 
quencies  of  interest.  Detonation  of  Composition  C-4  supplied  a  nonlinear, 
low-frequency-rich  source.  The  propane  cannon  provided  a  more  linear 
impulsive  source,  with  ample  low-frequency  energy.  The  loudspeaker  pro¬ 
vided  a  continuous-type  source  for  a  broad  range  of  frequencies.  The  pro¬ 
pane  cannon  and  loudspeaker  were  portable,  allowing  for  the  use  of  four 
separate  source  locations  with  minimal  fuss.  The  pistol  was  used  to  char¬ 
acterize  ground  impedance.  Limited  data  are  available  for  this  sound 
source. 

Scope 

This  report  is  intended  to  provide  a  full  description  of  the  data  set  from  the 
experiment  described  within  it.*  A  quality  evaluation  of  the  data  is  in¬ 
cluded  in  this  report,  as  well  as  a  description  of  data  reduction.  Analysis  of 
the  data  is  reserved  for  a  later  report. 


*  Data  sets  are  available  on  DVD  upon  request.  Send  requests  to  Dr.  Michelle  Swearingen,  PO  Box  9005, 
Champaign,  IL  61826-9005  or  email  her  at  Michelle.E.Swearingen@erdc.usace.army.mil. 
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2  Test  Overview 

The  impact  of  a  forest  edge  on  acoustic  propagation  is  currently  unknown. 
It  is  hypothesized  that  the  forest  edge  contributes  significant  attenuation 
when  in  a  propagation  path.  This  attenuation  will  be  caused  by  changes  in 
the  microclimate,  ground,  and  biomass  compared  to  an  open  area.  This 
basic  research  project  is  designed  to  isolate  and  study  each  factor  that 
could  influence  the  acoustic  propagation. 

The  experimental  measurements  were  conducted  on  Range  79  at  Camp 
Ripley,  MN.  Figure  1  contains  an  aerial  photograph  of  the  site  with  sensor 
locations  marked.  Figure  2  is  a  schematic  of  the  sensor  layout. 

Acoustical  data  recording  and  reduction 

Acoustical  data  were  recorded  on  two  model  DL750  16-channel  16-bit  Yo- 
kogawa  digital  recording  oscilloscopes,  using  type  701251  voltage  input 
cards.  The  four  four-element  microphone  arrays  were  recorded  on  Trailer 
scope  (located  in  the  trailer),  and  the  other  microphones  (S1-5,  L1-3)  were 
recorded  on  the  Forest  scope  (located  in  the  forest). 

Composition  C-4  records  were  5  seconds  long.  The  Trailer  scope  was  trig¬ 
gered  by  microphone  A2a;  the  Forest  scope  was  triggered  by  microphone 
S2.  In  each  case,  the  trigger  signal  contains  a  1  second  pre-trigger  and 
4  seconds  post-trigger.  The  sampling  rate  was  set  to  look  samples/second 
for  each  channel.  Other  microphone  channels  are  recorded  simultaneously 
and  show  the  blast  at  increasingly  later  times  with  increasing  distance 
from  the  blast  point. 

Propane  cannon  records  were  in  the  same  format  and  sampling  rate  as  the 
Composition  C-4  records.  However,  the  trigger  microphone  varied  accord¬ 
ing  to  blast  point.  Table  1  indicates  the  blast  point  and  trigger  microphone 
for  each  oscilloscope. 

Loudspeaker  records  were  manually  started  and  stopped  at  each  oscillo¬ 
scope  in  order  to  accommodate  the  lower  received  levels  and  to  avoid  re¬ 
cording  only  wind  noise.  Records  were  100  seconds  in  length  with  a  sam¬ 
pling  rate  of  50,000  samples/second. 
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Pistol  shot  data  were  only  recorded  on  four  1/2-in.  microphones,  placed  on 
the  ground  at  locations  A1-A4.  These  data  were  recorded  with  a  record 
length  of  2  seconds  and  a  sampling  rate  of  5,000  samples/second.  Trigger¬ 
ing  was  accomplished  manually. 


Ki 


Key: 

ropane  Cannon/C4 
Prop.  Can./Speaker 
Mic  Array 
WjJ  IO  Single  Mic 
•  /  □  Met  Tower 

'  > 


ZD  Meters 
150 


BP4 


Figure  1.  Aerial  view  of  test  site  with  sources,  sensors,  and  meteorological  instrumentation  locations  overlaid. 
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Figure  2.  Schematic  of  sensor  layout.  Green  hashmarked  box  represents  a  portion  of  the 
forest.  Forest  continues  in  the  +x  direction. 


Table  1.  Trigger  microphone  used  for 
each  propane  cannon  blast  point. 


Blast  Point 

Trailer 

Forest 

BP1 

Ala 

SI 

BP2 

Ala 

SI 

BP3 

A4a 

LI 

BP4 

A4a 

L3 
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Acoustical  data  in  all  cases  were  reduced  using  Matlab  scripts  (attached  in 
the  appendix).  As  the  data  were  all  recorded  digitally  as  integer-time  se¬ 
ries,  the  only  transformation  necessary  was  to  convert  the  files  from  the 
binary  oscilloscope  format  into  floating-point  Matlab  format  as  voltage¬ 
time  series.  Waveforms  remain  in  volts  but  are  converted  to  Pascals  by  ap¬ 
plying  pressure  voltage  sensitivities  within  plotting  scripts. 

List  of  signature  records 

The  inserted  Excel  spreadsheet  contains  the  following  information  for  all 
recordings:  Yokogawa  scope  file  number,  Date,  Time  (l  hour  fast),  Source 
Type,  Source  Location,  Array  Configuration,  Matlab  file  name/location, 
and  Comments.  The  comments  contain  information  on  data  quality.  Two 
worksheets  are  in  the  file,  one  for  the  Trailer  Scope  and  one  for  the  Forest 
Scope.  Table  2  lists  the  microphone  locations  and  appropriate  scope  and 
channel. 


Table  2.  Microphone  location,  recording 
oscilloscope,  and  channel  matrix. 


Microphone 

Scope  Loc. 

Channel 

Ala 

Trailer 

CHI 

Alb 

Trailer 

CH2 

Ale 

Trailer 

CH3 

Aid 

Trailer 

CH4 

A2a 

Trailer 

CH5 

A2b 

Trailer 

CH6 

A2c 

Trailer 

CH7 

A2d 

Trailer 

CH8 

A3a 

Trailer 

CH9 

A3b 

Trailer 

CH10 

A3c 

Trailer 

CH11 

A3d 

Trailer 

CH12 

A4a 

Trailer 

CH13 

A4b 

Trailer 

CH14 

A4c 

Trailer 

CH15 

A4d 

Trailer 

CH16 

LI 

Forest 

CHI 

L2 

Forest 

CH2 

L3 

Forest 

CH3 

SI 

Forest 

CH4 

S2 

Forest 

CH5 

S3 

Forest 

CH6 

S4 

Forest 

CH7 

S5 

Forest 

CH8 
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The  embedded  spreadsheet  (also  included  on  the  DVD)  contains  a  list  of 
all  recordings.  Quality  control  on  the  data  has  been  completed  in  the  sense 
that  all  missing  and  completely  unusable  signals  are  marked.  Noisy  signals 
are  not  necessarily  marked.  The  following  abbreviations  are  used: 

Pre-Cal:  Pre-calibration  signal 

Post-Cal:  Post-calibration  signal 

C-4:  Composition  C-4  explosive 

PC:  Propane  Cannon 


N:\Research\ 
Acoustics\Camp  Riple 


Specifications  on  sources 


C-4 

1  Brick 

1.25  lb  (0.57  kg) 

0.45  Pistol  Blank 

Overall  Length 

0.89  in.  (22.6  mm) 

Case  Length 

0.89  in.  (22.6  mm) 

Cartridge  Weight 

7-6  g 

Case  Weight 

5-8  g 

Powder  Weight 

31  gr 

Powder  Type 

Black,  FFFFG,  CCI 

Primer  Type 

Magnum 

Manufacturer 

Custom  made 

Propane  Cannon 

Manufacturer 

Reed-Joseph  International  Co. 

Model 

SCARE  AWAY  M-4  Cannon 

Loudspeaker 

Manufacturer 

Cerwin-Vega! 

Model 

PSX-122 

Amplifier 

Crown  Macro-Tech  600 

Configuration  details 

Tables  3  through  11  contain  sensor  and  source  locations,  as  well  as  sensor 
details  and  changes/issues. 
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Table  3.  Locations  of  microphones  and  blast  points.  Distances  are  to  BP1. 


Location  Name 

Coordinates 

Distance  to  BP1  (m) 

BP1 

(0,0) 

0 

BP2 

(-88.85,0) 

88.85 

BP3 

(0,175) 

175 

BP4 

(0,313) 

313 

Al(a-d) 

(0,25) 

25 

A2(a-d) 

(0,50) 

50 

A3(a-d) 

(0,75) 

75 

A4(a-d) 

(0,100) 

100 

SI 

(-10.6,22.7) 

25 

S2 

(-23.2,50) 

55.1 

S3 

(-23.2,75) 

78.5 

S4 

(-84.5,196) 

213.4 

S5 

(-97,196) 

218.7 

LI 

(0,150) 

150 

L2 

(0,263) 

263 

L3 

(0,288) 

288 

Table  4.  Initial  set  of  meteorological  sensors  on  the  open  field  tower. 


Location 

Height  (m) 

Height  (ft) 

Sensor  Type 

Serial  # 

Open 

Ground 

Ground 

Soil  Moisture 

1001458 

Open 

N/A 

N/A 

Barometric  Pressure 

718072 

Open 

~1 

3.28 

Net  Radiometer 

730735 

Open 

2.59 

8.5 

Anemometer 

716311 

Open 

2.59 

8.5 

8-bit  temp/humidity 

725373 

Open 

2.59 

8.5 

12-bit  temperature 

732365 

Open 

5.18 

17 

Anemometer 

981368 

Open 

5.18 

17 

8-bit  temp/humidity 

725365 

Open 

5.18 

17 

12-bit  temperature 

732762 

Open 

7.77 

25.5 

Anemometer 

991685 

Open 

7.77 

25.5 

8-bit  temp/humidity 

725370 

Open 

7.77 

25.5 

12-bit  temperature 

1000055 

Open 

10.36 

34 

Anemometer 

713662 

Open 

10.36 

34 

8-bit  temp/humidity 

725372 

Open 

10.36 

34 

12-bit  temperature 

732773 

Open 

12.95 

42.5 

Anemometer 

991684 

Open 

12.95 

42.5 

8-bit  temp/humidity 

725363 
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Location 

Height  (m) 

Height  (ft) 

Sensor  Type 

Serial  # 

Open 

12.95 

42.5 

12-bit  temperature 

732766 

Open 

~1 

3.28 

Rain  Gauge 

668179 

Open 

N/A 

N/A 

Data  Logger  -temperature 

724611 

Open 

N/A 

N/A 

Data  Logger  -  wind 

724610 

Open 

N/A 

N/A 

Data  Logger  -  other 

724615 

Table  5.  Initial  set  of  meteorological  sensors  on  the  forest  edge  tower. 


Location 

Height  (m) 

Height  (ft) 

Sensor  Type 

Serial  # 

Edge 

Ground 

Ground 

Soil  Moisture 

1001459 

Edge 

N/A 

N/A 

Barometric  Pressure 

718077 

Edge 

~1 

3.28 

Net  Radiometer 

730732 

Edge 

2.59 

8.5 

Anemometer 

981384 

Edge 

2.59 

8.5 

8-bit  temp/humidity 

725371 

Edge 

2.59 

8.5 

12-bit  temperature 

732761 

Edge 

5.18 

17 

Anemometer 

981383 

Edge 

5.18 

17 

8-bit  temp/humidity 

725356 

Edge 

5.18 

17 

12-bit  temperature 

732759 

Edge 

7.77 

25.5 

Anemometer 

713665 

Edge 

7.77 

25.5 

8-bit  temp/humidity 

725360 

Edge 

7.77 

25.5 

12-bit  temperature 

732760 

Edge 

10.36 

34 

Anemometer 

981385 

Edge 

10.36 

34 

8-bit  temp/humidity 

999197 

Edge 

10.36 

34 

12-bit  temperature 

732764 

Edge 

12.95 

42.5 

Anemometer 

713660 

Edge 

12.95 

42.5 

8-bit  temp/humidity 

725367 

Edge 

12.95 

42.5 

12-bit  temperature 

732768 

Edge 

~1 

3.28 

Rain  Gauge 

668173 

Edge 

N/A 

N/A 

Data  Logger  -temperature 

535230 

Edge 

N/A 

N/A 

Data  Logger  -  wind 

550714 

Edge 

N/A 

N/A 

Data  Logger  -  other 

595411 

Table  6.  Initial  set  of  meteorological  sensors  on  the  forest  tower. 


Location 

Height  (m) 

Height  (ft) 

Sensor  Type 

Serial  # 

Forest 

Ground 

Ground 

Soil  Moisture 

986131 

Forest 

N/A 

N/A 

Barometric  Pressure 

718071 

Forest 

~1 

3.28 

Net  Radiometer 

730731 

Forest 

2.59 

8.5 

Anemometer 

713657 
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Location 

Height  (m) 

Height  (ft) 

Sensor  Type 

Serial  # 

Forest 

2.59 

8.5 

8-bit  temp/humidity 

725358 

Forest 

2.59 

8.5 

12-bit  temperature 

732763 

Forest 

5.18 

17 

Anemometer 

713658 

Forest 

5.18 

17 

8-bit  temp/humidity 

725357 

Forest 

5.18 

17 

12-bit  temperature 

732767 

Forest 

7.77 

25.5 

Anemometer 

965682 

Forest 

7.77 

25.5 

8-bit  temp/humidity 

725359 

Forest 

7.77 

25.5 

12-bit  temperature 

732770 

Forest 

10.36 

34 

Anemometer 

713666 

Forest 

10.36 

34 

8-bit  temp/humidity 

725366 

Forest 

10.36 

34 

12-bit  temperature 

732772 

Forest 

12.95 

42.5 

Anemometer 

716310 

Forest 

12.95 

42.5 

8-bit  temp/humidity 

725364 

Forest 

12.95 

42.5 

12-bit  temperature 

732769 

Forest 

~1 

3.28 

Rain  Gauge 

668171 

Forest 

N/A 

N/A 

Data  Logger  -temperature 

724627 

Forest 

N/A 

N/A 

Data  Logger  -  wind 

724625 

Forest 

N/A 

N/A 

Data  Logger  -  other 

724608 

Table  7.  Issues  with  meteorological  sensors  during  the  experiment. 


Location 

Height 

(m) 

Sensor  Type 

Serial  # 

Date 

Failed 

Time 

Failed 

Replaced? 

Open 

N/A 

Data  Logger 

724611 

unknown 

unknown 

724619  on  06/24/06 
15:55 

Edge 

2.59 

8-bit 

temp/humidity 

725371 

N/A 

N/A 

No.  Manufacturer 
evaluation  indicates 
failure 

Forest 

12.95 

12-bit 

temperature 

732769 

24  June 

09:20 

No.  Manufacturer 
evaluation  indicates  no 
defect 

Open 

12.95 

Anemometer 

991684 

24  June 

05:10 

No.  Manufacturer 
evaluation  indicates  no 
defect. 
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Table  8.  Composition  C-4  microphone  list.  Distances  are  from  BP1. 


Location 

Distance  (m) 

Mic.  Type 

Microphone  # 

Power  Supply  # 

Changed? 

Ala 

25 

Blast  Pencil 

5872 

8014 

N 

Alb 

25 

Blast  Pencil 

5870 

8012 

N 

Ale 

25 

Blast  Pencil 

5875 

8197 

N 

Aid 

25 

Blast  Pencil 

5869 

8013 

N 

A2a 

50 

1/8  in. 

2515896 

1466052 

N 

A2b 

50 

1/8  in. 

45967 

1466052 

N 

A2c 

50 

1/8  in. 

56385 

1466063 

N 

A2d 

50 

1/8  in. 

56396 

1466063 

N 

A3  a 

75 

1/8  in. 

2515907 

1533630 

N 

A3b 

75 

1/8  in. 

2515897 

1533630 

N 

A3c 

75 

1/8  in. 

2515914 

1466060 

N 

A3d 

75 

1/8  in. 

2515901 

1466060 

N 

A4a 

100 

1/4  in. 

991054 

1533663 

N 

A4b 

100 

1/4  in. 

38681 

1533663 

N 

A4c 

100 

1/4  in. 

41576 

2424025 

N 

A4d 

100 

1/4  in. 

991206 

2424025 

N 

SI 

25 

Blast  Pencil 

5871 

8015 

N 

S2 

55.1 

1/8  in. 

56378 

555809 

N 

S3 

78.5 

1/8  in. 

2515888 

1466091 

N 

S4 

213.4 

1/4  in. 

38671 

1533661 

N 

S5 

218.7 

1/4  in. 

35968 

1533632 

N 

LI 

150 

1/4  in. 

41574 

1466051 

N 

L2 

263 

1/4  in. 

35976 

936299 

N 

L3 

288 

1/4  in. 

35979 

1466048 

N 

Table  9.  Propane  cannon  microphone  list.  Distances  are  from  BP1. 


Location 

Distance  (m) 

Mic.  Type  (in.) 

Microphone  # 

Power  Supply  # 

Changed? 

Ala 

25 

1/4 

41584 

761782 

N 

Alb 

25 

1/4 

41572 

761782 

N 

Ale 

25 

1/4 

41571 

1533633 

N 

Aid 

25 

1/4 

35975 

1533633 

N 

A2a 

50 

1/4 

41577 

1466052 

N 

A2b 

50 

1/4 

318293 

1466052 

N 
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Location 

Distance  (m) 

Mic.  Type  (in.) 

Microphone  # 

Power  Supply  # 

Changed? 

A2c 

50 

1/4 

48592 

1466063 

N 

A2d 

50 

1/4 

41579 

1466063 

N 

A3  a 

75 

1/4 

41520 

1533630 

N 

A3b 

75 

1/4 

38674 

1533630 

N 

A3c 

75 

1/4 

35966 

1466060 

N 

A3d 

75 

1/4 

52399 

1466060 

Y 

A4a 

100 

1/4 

991054 

1533663 

N 

A4b 

100 

1/4 

38681 

1533663 

Y 

A4c 

100 

1/4 

41576 

2424025 

N 

A4d 

100 

1/4 

991206 

2424025 

N 

SI 

25 

1/4 

52398 

1466055 

N 

S2 

55.1 

1/4 

41578 

555809 

N 

S3 

78.5 

1/4 

41573 

1466091 

N 

S4 

213.4 

1/4 

38671 

1533661 

N 

S5 

218.7 

1/4 

35968 

1533632 

N 

LI 

150 

1/4 

41574 

1466051 

N 

L2 

263 

1/4 

35976 

936299 

N 

L3 

288 

1/4 

35979 

1466048 

N 

Table  10.  Loudspeaker  microphone  list.  Distances  are  from  BP1. 


Location 

Distance  (m) 

Mic.  Type  (in.) 

Microphone  # 

Power  Supply  # 

Changed? 

Ala 

25 

1/4 

41584 

761782 

N 

Alb 

25 

1/4 

41572 

761782 

N 

Ale 

25 

1/4 

41571 

1533633 

N 

Aid 

25 

1/4 

35975 

1533633 

N 

A2a 

50 

1/2 

1857571 

1466052 

N 

A2b 

50 

1/2 

1357654 

1466052 

N 

A2c 

50 

1/2 

1881050 

1466063 

N 

A2d 

50 

1/2 

2335415 

1466063 

N 

A3  a 

75 

1/2 

1881043 

1533630 

N 

A3b 

75 

1/2 

1881046 

1533630 

N 

A3c 

75 

1/2 

1881049 

1466060 

N 

A3d 

75 

1/2 

2068999 

1466060 

N 

A4a 

100 

1/2 

1857589 

1533663 

N 
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Location 

Distance  (m) 

Mic.  Type  (in.) 

Microphone  # 

Power  Supply  # 

Changed? 

A4b 

100 

1/2 

1848766 

1533663 

N 

A4c 

100 

1/2 

1848782 

2424025 

N 

A4d 

100 

1/2 

1783511 

2424025 

N 

SI 

25 

1/2 

52398 

1466055 

N 

S2 

55.1 

1/2 

1838395 

555809 

N 

S3 

78.5 

1/4 

41573 

1466091 

N 

S4 

213.4 

1/4 

38671 

1533661 

N 

S5 

218.7 

1/4 

35968 

1533632 

N 

LI 

150 

1/4 

41574 

1466051 

N 

L2 

263 

1/4 

35976 

936299 

N 

L3 

288 

1/4 

35979 

1466048 

N 

Table  11.  Issues  with  microphones  during  propane  cannon  test. 


Location 

Mic  Type  (in.) 

Serial  # 

Date  Failed 

Time  Failed 

Replaced? 

A3d 

1/4 

52399 

24  June 

13:16 

2239253  on  6/24/06  14:00 

A4b 

1/4 

38681 

24  June 

13:16 

2239254  on  6/24/06  14:00 

Microphone  array 

The  microphones  set  up  at  locations  A1-A4  all  included  a  four- 
microphone  array  set  in  a  minimum  baseline  redundancy  array  configura¬ 
tion.  A  schematic  of  this  configuration  is  shown  in  Figure  3.  Figure  4 
shows  a  microphone  array  in  a  parallel  configuration. 


Figure  3.  Minimum  baseline-redundancy  array.  This  array  provides  six 
baselines  per  four  microphones.  Locations  2,  3,  and  5  are  not  used. 
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Figure  4.  A  microphone  array  is  shown  in  the  parallel  configuration. 


In  this  configuration,  the  following  distances  were  used:  (a)  66.5  cm  to  the 
left  of  the  center  point,  (b)  41  cm  to  the  left  of  the  center  point,  (c)  34.5  cm 
to  the  right  of  the  center  point,  and  (d)  84  cm  to  the  right  of  the  center 
point.  In  the  perpendicular  configuration,  (a)  was  to  the  left  of  the  propa¬ 
gation  line.  In  the  parallel  configuration,  (a)  was  closest  to  the  source 
point.  In  the  vertical  configuration,  (a)  was  the  highest  microphone.  The 
height  was  set  so  that  in  the  perpendicular  and  parallel  configurations,  all 
microphones  were  1.4  m  above  the  ground  surface. 

Microphone  issues 

During  the  experiment,  some  of  the  microphones  failed  to  produce  a  clean 
calibration  signal.  When  this  problem  occurred,  the  microphone,  pre¬ 
amplifier,  and/or  power  supply  were  changed  and  the  new  serial  number 
of  each  new  component  was  used  for  all  of  the  following  recordings  (see 
Table  11).  In  no  instances  did  a  microphone  fail  during  a  test.  All  issues 
were  resolved  during  pre-calibrations,  and  none  were  found  during  post¬ 
calibrations. 

Meteorological  sensor  issues 

During  the  experiment,  the  temperature  data  logger  on  the  open  field 
tower  failed.  It  was  replaced  on  24  June  and  all  temperature  and  tempera¬ 
ture/humidity  sensors  checked  in  as  working.  The  failed  data  logger  was 
sent  to  the  manufacturer  in  the  hopes  that  the  data  could  be  recovered. 
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However,  the  manufacturer  was  not  able  to  find  anything  wrong  with  the 
logger  or  retrieve  any  data.  This  particular  data  logger  was  facing  the 
Composition  C-4  blasts  from  a  distance  of  approximately  50  m.  The  prox¬ 
imity  to  the  blast  is  a  possible  cause  for  the  failure. 

During  the  experiment,  there  were  occasional  dropped  recordings.  It  is 
possible  that  the  l-minute  polling  delay  on  the  sensors  was  not  long 
enough  for  the  number  of  sensors  in  the  data  loggers.  Two  of  the  sensors 
(listed  in  Table  7)  failed  during  the  experiment.  They  were  sent  to  the 
manufacturer  and  found  to  be  free  of  defects.  The  reason  for  failure  is  un¬ 
known.  One  additional  sensor  (also  listed  in  Table  7)  appeared  to  give 
questionable  data.  It  was  sent  to  the  manufacturer  and  determined  to  be 
faulty  and  not  repairable. 

Sensor  calibration 

Microphones  were  calibrated  both  before  and  after  each  recording  session 
using  a  B&K  type  4228  pistonphone  calibrator  with  appropriate  adapters. 
The  4228  calibrators  produce  an  equivalent  level  equal  to  124  dB  root 
mean  square  (peak  sound  pressure  level  equal  to  127  dB)  at  250  Hz.  The 
calibration  procedure  required  that  the  calibrator  and  microphone  lie  on 
the  ground  to  minimize  extraneous  vibrations.  Clean  5-second  calibration 
signals  at  a  sampling  rate  of  100,000  samples/second  were  recorded  for 
each  microphone. 

The  calibration  procedure  for  the  blast  gauges  was  much  more  primitive. 
Blast  gauges  were  struck  with  a  hard  object  to  induce  a  signal.  If  a  clear 
signal  was  detected,  the  blast  gauge  was  considered  operational.  Calibra¬ 
tion  values  for  each  blast  gauge  were  obtained  from  calibration  sheets  pro¬ 
vided  during  the  last  manufacturer  calibration.  Previous  experience  indi¬ 
cates  that  the  calibration  values  of  the  blast  gauges  do  not  change 
significantly  over  time. 

Representative  acoustic  data 

Figures  5  through  8  show  samples  of  the  pressure-time  series  data  col¬ 
lected  during  this  experiment.  In  each  figure,  the  source  location  is  BPi. 
Microphone  locations  are  listed  on  the  figures.  All  waveforms  are  cali¬ 
brated  pressure  signatures.  Each  of  the  four  source  types  is  represented: 
Composition  C-4,  Propane  Cannon,  Loudspeaker,  and  Pistol  Shot. 


C4  sample  from  2006/06*22  at  13:02:19.59 
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Figure  5.  Sample  pressure-time  histories  from  a  single  shot  of  Composition  C-4  at  shot  location  BP1. 
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Propane  Cannon  sample  from  2006/06/23  at  12:50:31.33 
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Figure  6.  Sample  pressure-time  histories  from  a  single  propane  cannon  shot  at  shot  location  BP1. 
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Loudspeaker  sample  from  2006/06/25  at  14:26:20.21. 
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Figure  7.  Sample  pressure-time  histories  from  a  single  loudspeaker  run  at  shot  location  BP1. 
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Pistol  sample  from  2006/06/22  at  16:15:03.28. 
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Figure  8.  Sample  pressure-time  histories  from  a  single  pistol  shot. 
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Representative  meteorological  data 

The  following  three  tables  show  samples  of  the  meteorological  data  col¬ 
lected  during  this  experiment.  The  three  tower  locations  are  open,  edge, 
and  forest.  Data  are  presented  for  temperature  and  wind  speed  for  a  single 
time.  Relative  humidity  and  wind  direction  were  also  recorded.  During  the 
experiment,  data  from  each  sensor  were  recorded  once  every  minute. 
Every  5  minutes  the  data  logger  performed  an  average  over  the  previous 
five  recordings  and  saved  the  5-minute  average  to  its  storage  device.  All  of 
the  data  available  from  this  experiment  are  in  the  form  of  5-minute  aver¬ 
ages. 


Open  Edge  Forest 


Figure  9.  Temperature  profiles  at  all  three  towers  at  the  same  time. 
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Open  Edge  Forest 


Figure  10.  Wind  speed  profiles  from  all  three  towers  at  the  same  time. 

Forest/plantation  characterization 

Characterization  of  the  forest  vegetation  will  be  conducted  at  two  levels. 
The  primary  level  is  the  forest  edge  (Edge  Characterization);  the  secondary 
is  the  forest  properties  in  general  (Stand  Characterization).  Edge  charac¬ 
terization  will  include  the  measurement  and  characterization  of  all  woody 
and  foliage  components  from  the  edge  of  the  plantation  (where  the  trees 
meet  the  adjacent  grassland)  to  the  where  the  canopy  height/density  nor¬ 
malizes  (non-edge  plantation).  There  is  a  lack  of  woody/ shrub  component 
along  the  edge  so  there  is  no  need  to  concentrate  on  forest  under  story 
components.  Stand  characterization  will  mainly  focus  on  the  corridor  of 
trees  along  the  microphone  array,  with  secondary  measurements  taken  to 
account  for  the  areas  not  on  in  the  corridor.  The  forest  consisted  entirely 
of  Red  Pine  ( Pinus  resinosa).  Figures  11-14  show  varying  views  of  the  for¬ 
est. 


ERDC  SR-07-3 


22 


Figure  11.  Forest  edge  photographed  from  the  open  field. 


Figure  12.  Forest  edge  photographed  from  inside  the  forest. 
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Figure  13.  Typical  forest  floor  at  the  test  site. 


Figure  14:  Forest  interior  with  microphone  set  in  the  center  of  the  photograph. 
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Edge  characterization 

The  structure  of  the  plantation’s  border  between  the  trees  and  adjacent 
grassland  is  very  simple.  There  are  no  shrubs/woody  species  along  the 
edge  other  than  trees,  there  is  no  tree  recruitment  (all  trees  are  the  same 
age),  and  there  are  minimal  numbers  of  grasses  and  nonwoody  vegetation. 
This  structure  is  mainly  due  to  the  active  management  practices  (fire)  used 
to  manage  the  adjacent  grassland.  Additionally,  all  trees  are  in  distinct 
rows. 

1.  Individual  tree  measurements: 

Individual  trees  in  the  edge  were  measured  according  to  their  grouping  in 
rows  parallel  from  the  edge  (i.e.,  first  row  measured,  second  row,  etc.).  In¬ 
formation  recorded: 

a.  Diameter  at  Breast  Height  (DBH) 

b.  Species 

c.  Location 

d.  Height 

e.  Height  from  ground  to  canopy 

2.  Tree  component  measurements: 

For  each  row  of  trees,  the  region  from  the  edge  to  a  point  where  the  canopy 
normalizes  the  structural  components  of  the  trees  was  quantified  (volume 
of  foliage,  branches,  and  twigs  (fine  branches).  A  representative  sample  of 
trees  was  measured  to  establish  parameters  for  the  density  and  volumes  of 
these  materials. 

Selected  trees  were  chosen,  their  limbs  were  removed  up  to  the  canopy 
height,  and  materials  were  separated  into  branches,  foliage,  and  twigs. 
Branch  diameters  were  measured  at  both  ends  and  at  the  center  with  a 
caliper.  Additionally,  heights  were  measured  (up  to  the  maximum  height 
of  limb  removal)  and  a  measurement  was  taken  to  account  for  the  total  ra¬ 
dius  of  the  on-tree  foliage.  Taking  these  steps  allowed  for  determining  the 
area  of  foliage.  The  three  components  separated  from  the  tree  were  meas¬ 
ured  to  account  for  volume  of  foliage,  branches,  and  twigs.  These  volumes 
were  determined  by  bulk  immersion  in  a  30-gallon  trash-can  filled  with 
water  and  measuring  the  displaced  volume. 
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3 .  Digital  imagery  analysis : 

In  addition  to  the  physical  measurements  of  forest  edge  components;  the 
density  of  vegetation  on  the  edge  was  measured  using  digital  imagery 
analysis.  A  white  cloth  sheet  was  placed  at  intervals  behind  trees  along  the 
edge  and  a  photograph  was  taken.  These  intervals  were  selected  based  on 
the  parallel  rows  of  trees  along  the  edge.  Analysis  software  will  be  used  to 
delineate  the  percentage  of  foliage  and  woody  materials  visible  at  each  in¬ 
terval. 

Stand  characterization  measurements 

In  addition  to  characterizing  the  forest  edge,  sampling  was  conducted  to 
characterize  stand  parameters.  A 100%  survey  was  conducted  along  the 
microphone/instrument  array.  Based  on  estimates  of  the  time  interval 
needed  to  account  for  sound  reverberations  from  C-4  charges,  a  corridor 
was  delineated  along  the  array  and  all  trees  were  tallied.  The  following  in¬ 
formation  was  recorded  for  all  trees  within  the  corridor: 

1.  Diameter  at  Breast  Height  (DBH)  (DBH  Tape) 

2.  Species 

3.  Location 

(Transect  Width:  20-m  width  [10  m  each  side  of  transect  line].) 

Stand  parameters  within  the  corridor  were  measured  to  generally  charac¬ 
terize  the  forested  environment.  These  parameters  (and  methods  of  meas¬ 
urement)  are: 

1.  Crown  Closure  (optical  point  sampler) 

2.  Tree  Height  (clinometer) 

3.  Canopy  Height  (clinometer) 

4.  Basal  Area  (calculation) 

5.  Trees  per  Acre  (calculation). 

The  forested  area  outside  the  measurement  corridor  was  sampled  using  a 
systemic  random  methodology.  Variable  radius  plots  were  established 
within  these  areas  using  a  10  basal  area  factor  prism.  The  same  individual 
and  stand  level  parameters  were  measured  within  these  plots. 

Additionally,  the  volume  density  of  dead  branches  still  attached  to  trunks 
was  estimated  in  the  plantation  using  methods  similar  to  those  for  compo- 


ERDC  SR-07-3 


26 


nent  measurement  on  the  edge.  Most  of  the  trees  within  the  plantation 
have  retained  their  dead  understory  branches. 

Forest  environment 

The  forest  was  characterized  thoroughly.  The  following  section  contains 
details  of  these  measurements. 

General  forest  stats 

Basal  Area:  125  sq  ft/acre 

Trees  Per  Acre  (TPA):  327 

Average  DBH:  8.2  in. 

Array  transect  stats 

Transect  was  20  m  by  125  m  centered  on  centerline  of  array.  The  starting 
(o  meter  mark)  point  was  the  microphone  stand  on  the  forest  edge. 

Average  DBH:  7.7  in. 

TPA  (along  total  array):  429 

Tree  Heights:  36  ft 
Height  to  Canopy:  19  ft 

Tree  Height  EDGE:  26  ft 
Canopy  Height  EDGE:  5  ft 

Litter  characteristics 

Litter  was  sampled  at  each  of  the  microphone  points  along  the  array  up  to 
Li.  Additional  samples  were  taken  at  midpoints  between  the  microphones. 
One-meter  square  areas  of  litter  were  collected  to  measure  volume. 
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Table  12.  Litter  layer  characteristics. 


Sample  Point 

Depth  (mm) 

Volume  (cubic  cm) 

A2 

13 

946 

31 

2,120 

A3 

30 

5,413 

41 

6,852 

A4 

36 

5,186 

17 

3,066 

LI 

31 

4,732 

Woods  Blast  Point 

28 

2,839 

Crown  closure 

Crown  closure  was  measured  at  25-m  intervals  starting  25  m  into  the 
woods  from  A2.  An  optical  device  was  used. 


Table  13.  Crown  closure. 


Point 

Closure  (%) 

25  meters  in 

55 

A3 

70 

60 

A4 

70 

50 

flag 

40 

70 

LI 

75 

65 

Blast  Point  Woods 

80 

Tree  biomass  data 

Several  trees  were  destructively  sampled  to  determine  their  characteris¬ 
tics.  The  following  information  describes  these  characteristics. 

Branches  removed  up  to  10  ft  from  trees  within  the  woods;  15  ft  for  edge 
tree. 
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Tree  1 

Branch  radius:  7.5  ft 
DBH:  8.9  in. 

Ht:  41  ft 
Canopy:  20  ft 

Dead  Branches:  Number:  22 

Avg.  Length:  7.03  ft 

Avg.  Fat  Diameter:  1.29  in. 

Avg.  Mid-point  Diameter:  0.83  in. 

Avg.  End-point  Diameter:  0.38  in. 

Biomass  Volumes:  Branches:  10,410  cubic  centimeter  [cc]) 
Twigs:  11,356  cc 


Tree  2 

Branch  radius:  4.0  ft 
DBH:  6.3  in. 

Ht:  33  ft 
Canopy:  16  ft 

Dead  Branches:  Number:  25 

Avg.  Length:  5.12  ft 

Avg.  Fat  Diameter:  79  in. 

Avg.  Mid-point  Diameter:  0.48  in. 
Avg.  End-point  Diameter:  0.17  in. 

Biomass  Volumes:  Branches:  5,678  cc 
Twigs:  2,839  cc 


Tree  3 

Branch  radius:  6.0  ft 
DBH:  8.7  in. 

Ht:  45  ft 
Canopy:  26  ft 

Dead  Branches:  Number:  31 

Avg.  Length:  5.08  ft 

Avg.  Fat  Diameter:  1.12  in. 
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Avg.  Mid-point  Diameter:  0.74  in. 
Avg.  End-point  Diameter:  0.43  in. 

Biomass  Volumes:  Branches:  11,356  cc 
Twigs:  8,517  cc 


Tree  4 

Branch  radius:  9.0  ft 
DBH:  9.4  in. 

Ht:  38  ft 
Canopy:  20  ft 

Dead  Branches:  Number:  39 

Avg.  Length:  5.07  ft 

Avg.  Fat  Diameter:  0.95  in. 

Avg.  Mid-point  Diameter:  0.65  in. 
Avg.  End-point  Diameter:  0.39  in. 

Biomass  Volumes:  Branches:  13,249  cc 
Twigs:  12,303  cc 

Edge  Tree 

Branch  radius:  7.3  ft  (-3.5  at  15  in.) 

DBH:  8.0  in. 

Ht:  26  ft 
Canopy:  5  ft 
Dead  Branches: 


Live  Branches: 


Number:  10 

Avg.  Length:  6.32  ft 

Avg.  Fat  Diameter:  1.3  in. 

Avg.  Mid-point  Diameter:  0.81  in. 
Avg.  End-point  Diameter:  0.38  in. 

Number:  21 

Avg.  Length:  7.86  ft 

Avg.  Fat  Diameter:  1.6  in. 

Avg.  Mid-point  Diameter:  1.03  in. 
Avg.  End-point  Diameter:  0.48  in. 
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Biomass  Volumes:  Live  Branches:  31,230  cc 
Live  Twigs:  16,088  cc 
Needles:  16,088  cc 
Dead  Branches:  7,571  cc 
Dead  Twigs:  1,893  cc 

Given  the  edge  extends  5  m  in  (material  collected  within  15  ft  of  ground,  10 
trees  within  edge),  frontage  20  m  (100  m2  area): 

0.007  m3  (7287  cc)  of  canopy  material  per  square  meter 
0.003  m3  (3123  cc)  of  live  branch  material  per  square  meter 
0.002  m3  (1609  cc)  of  needle/live  twig  material  per  square  meter 

Archival  data 

All  binary  acoustic  files  from  the  Yokogawa  oscilloscope,  along  with  Mat- 
lab  scripts  to  reduce  the  data  into  Matlab  format  with  physical  units,  are 
included  on  the  available  data  disk  (see  footnote,  p  2).  Raw  weather  data 
in  Microsoft  Excel  format  are  included.  An  archive  of  site  photographs  and 
a  readme  file  that  describes  file  structure  and  simple  directions  on  using 
the  Matlab  scripts  are  also  included. 
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3  Summary  and  Conclusions 

The  experiment  was  successful  in  obtaining  excellent  quality  acoustical 
data  illuminating  the  effect  of  a  forest  edge  on  acoustic  propagation.  No 
acoustic  equipment  (e.g.,  sensors,  power  supplies,  recording  devices) 
failed  during  data-taking  sessions.  The  vast  majority  of  acoustic  signals 
recorded  are  easily  distinguishable  above  the  background  noise,  and  only 
one  data  point  out  of  several  thousand  was  clipped. 

The  weather  was  good  during  the  bulk  of  the  experiment,  with  rain  causing 
only  minimal  delays.  Winds  were  generally  light,  probably  a  desirable  at¬ 
tribute  for  this  measurement.  While  there  are  some  minor  deficiencies  in 
the  meteorological  data,  they  are  not  insurmountable.  Analysis  of  the  data 
will  determine  whether  the  experiment  has  answered  the  question  of 
“what  is  the  acoustic  effect  of  a  forest  edge.” 
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Appendix:  Matlab  Programs 

The  first  set  of  scripts  shown  below  was  used  to  convert  the  raw  data  from 
the  Yokogawa  to  a  Matlab-friendly  format.  These  scripts  were  developed 
for  previous  projects  using  the  same  Yokogawa  digital  oscilloscopes.  The 
data  from  the  Yokogawa  included  a  file  for  each  waveform  and  a  header 
with  pertinent  information  such  as  sampling  rate,  shot  time,  etc.  The  func¬ 
tions  organizeWaveforms, parseHeader,  writeSensitivities, 
and  f  indSensitivity  were  used  to  extract  this  information  and  reor¬ 
ganize  it  for  easy  analysis  in  Matlab. 

The  second  set  of  scripts  (beginning  p  43)  was  written  specifically  for  this 
project.  The  function  plotonl  is  used  to  view  all  microphone  signals  for  a 
single  shot  simultaneously.  This  function  is  useful  in  observing  the  change 
in  the  waveform  as  it  travels  through  the  air.  The  function  plotsingle 
plots  the  waveform  for  a  single  microphone.  This  function  offers  a  quick 
and  easy  way  to  view  any  one  channel’s  wave.  WaveandSpect  ra  similarly 
plots  a  single  waveform,  but  also  includes  a  one-third  octave  band  fre¬ 
quency  analysis.  The  functions  getPeaks,  and  plotPeaks  were  designed 
to  analyze  changes  in  peak  levels  versus  distance  from  the  source.  get¬ 
Peaks  determined  the  peak  levels  for  each  waveform  and  saved  them  in  a 
matrix  called  peaksT  able,  mat,  with  columns  indicating  which  source, 
blast  point,  and  array  configuration  were  used.  The  function  plotPeaks 
loads  the  matrix  created  by  getPeaks  and  then  creates  graphs  relating 
these  peak  levels  to  each  microphone’s  distance  from  the  source. 
PeaksPerShot  provides  a  view  of  the  peak  level  of  each  microphone  loca¬ 
tion  for  individual  shots  of  a  given  source,  source  location,  and  array  con¬ 
figuration. 

To  observe  effects  on  frequency  content  as  the  waves  traveled  through  the 
forest,  another  set  of  scripts  was  developed.  The  function  SELtest2  takes 
a  single  waveform  and  performs  a  one-third  octave  band  frequency  analy¬ 
sis.  SELavg2  provides  a  way  to  calculate  sound  exposure  levels  for  all 
shots  for  a  particular  source  type  and  blast  point.  The  results  can  either  be 
returned  for  each  shot  or  can  be  averaged  then  returned.  collectSEL2 
was  written  to  assemble  the  matrix  tobSEL_all .  mat  which  contains 
data  for  all  the  C-4  and  propane  cannon  shots,  along  with  their  overall  SEL 
and  one-third  octave  band  spectra.  SELperShot  works  like  PeaksPer¬ 
Shot,  except  that  it  displays  overall  SEL  values  instead  of  peak  pressures. 
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The  comments  accompanying  each  function  provide  more  detailed  expla¬ 
nations. 


function  organizeWaveforms (rootDir) 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Organizes  and  loads  the  .wvf  and  .hdr  files  corresponding  to  a  day  of 
%events  into  separate  directories,  saving  the  waveform  and  header 
%information  into  .mat  files 

%Written  for  experiments  in  Edgewood,  MD  June  2004  and  in  Blossom  Point,  MD 
%August  2004. 

%Version  1.0  completed  24  JUN  2004  by: 

%Tim  Eggerding 

% (t-eggerding@cecer.army.mil) 


%Update  history: 
%  12  AUG  2004: 


Use  of  fread  on  .wvf  file  changed  such  that  the  header 
data  is  read  initially  and  then  thrown  out,  and  the 
waveform  data  is  read  as  16  bit  integers  and  saved  as 
an  array  of  16  bit  integers.  fread  defaults  to  a 
return  type  of  a  double  array,  which  resulted  in  memory 
errors.  Casting  to  double  now  occurs  during  the 
process  of  breaking  the  .wvf  file  into  seperate 
waveforms  for  each  channel  and  converting  from  16  bit 
integers  to  64  bit  double  representing  accurate 
voltages . 

NOTE:  MATLAB  seems  to  run  out  of  memory  at  about 

50000000  doubles  in  one  array  due  to  contigious  memory 
issues.  If  more  than  50000000  samples  exist  in  a 
waveform,  more  than  one  array  must  be  created  in  order 
to  store  the  samples  with  double  precision. 


%  17  AUG  2004:  Changed  handling  of  filenames  to  allow  for  arbitrary 

%  length  file  names  as  opposed  to  restricting  to  Yokogawa-type 

%  4-character  file  names. 

%  24  SEP  2004:  Changed  the  variable  date  to  shotDate  to  avoid  conflicts 

%  built-in  MATLAB  function  date ( ) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 


%Ref erring  functions:  yokogawaDataGUI .m 


%Required  functions: 
%global  variables: 
%Inputs:  rootDir 


parseHeader .m 


Root  directory  for  the  day 


%Outputs : 


shotTimes .mat 
channels .mat 


Saves  an  array  of  times  of  events 

Saves  a  matrix  of  valid  channels  for  each  event 


%Creates  a  vector  listing  all  the  files  in  the  directory 
%numFiles  =  the  number  of  files  in  the  directory 


cd (rootDir) ; 

D  =  dir (rootDir) ; 
numFiles  =  length (D); 
shotTimeArray  =  [ ] ; 
channelArray  =  []; 

for  i  =  1: numFiles  %Look  through  all  the  files 

wvf Name  =  D(i) .name;  %Get  the  current  file  name 
if (~isempty (strfind (wvfName, ' WVF' ) ) )  %If  the  extension  is  .WVF 
cd (rootDir) 

nameLength=length (wvfName) -4; 
hdrName  =  [wvfName ( 1 : nameLength) 

.HDR  %f ile 

hdrID  =  f open (hdrName,  '  r'); 
if  hdrID  ==  -1 
turns  -1 

warning ([' Header  file  '  hdrName  '  for  waveform  file  '  wvfName  '  does  not  ex 

ist . ' ] ) ; 

else 

%Load  the  header  information  using  parseHeader .m 
[endian,  dataFormat,  dataOffset,  traceTotalNumber ,  traceNames,  blockSize, 
vResolution,  vOffset,  vDataType,  period,  dateCell,  timeCell]  =  parseHeader (hdrName) ; 


%take  the  filename  without  the  .WVF 

'.HDR'];  %Compose  the  name  for  the  corresponding 

%Open  the  .HDR  file 

%Warn  that  the  header  was  not  found  if  fopen  re- 
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if  strcmp (endian,  'Big') 

wvfID  =  f open (wvf Name, ' r ' ,  'ieee-be');  %Use  big  endian  encoding  if  nec¬ 
essary 

else 

wvfID  =  fopen  (wvfName,  'r',  'ieee-le');  %Use  little  endian  encoding  if 

necessary 

end 

shotTimeArray  =  [shotTimeArray;  timeCell ( 1 ) ] ;  %Append  the  time  of  the  current 
event  to  %the  shotTimeArray 

if (length (traceNames) <16) 
tNamel6  =  cell (1,16); 

tNamel6 (1 : length (traceNames) )  =  traceNames; 

else 


tNamel6  =  traceNames; 

end 

channelArray  =  [channelArray;  tNamel6] ; 
event  to  %the  channelArray 

DataType  =  char (vDataType (1) ) ; 
array  % (string) 

time  =  char (timeCell ( 1 ))  ; 

(string) 

shotDate  =  char (dateCell ( 1 ) )  ; 
array  % (string) 

dataSize  =  8*str2num (DataType (3)  )  ; 
dataString  =  num2str (dataSize) ; 
shotTime  =  [time (1:2)  time (4: 5) 

tion  to  %make  a  valid  directory  name 

use  a  %colon  in  directory  names) . 

if  ~ (exist (shotTime,  'dir')  ==  7) 
ready  been  %done 

shotTime 

mkdir (shotTime) ; 

event. 


%Append  the  channels  of  the  current 

%Put  the  data  type  into  a  character 

%Put  the  time  into  a  character  array 

%Put  the  date  into  a  character 

%Find  the  size  of  the  data  in  bits 
%Put  that  size  into  a  a  string 
time (7: 9)];  %Use  the  time  informa- 

%In  this  case,  hh.mm.ss  (you  can't 

%Make  sure  that  this  event  hasn't  al- 


%Make  the  subdirectory  to  hold  this 


cd (shotTime) ;  %and  enter  it. 

pack; 

headerData  =  f read (wvf ID, dataOff set/2 ,  ['*int'  dataString]);  %Read  the 

header  %data 


clear  headerData 

for  m  =  1 : length (traceNames)  %For  each  trace  in  the  header  file 

waveform  =  f read (wvfID,  blockSize (m) ,  ['*int'  dataString]); 

waveform  =  double (waveform) *vResolution (m) ; 
waveform  =  waveform  +  vOffset(l); 

dirName  =  char (traceNames (m) ) ;  %Set  the  directory  name  for  the  cur¬ 

rent  %channel's  data 

mkdir (dirName) ;  %Make  that  directory 

cd (dirName);  %Enter  that  directory 

hdrname  =sprintf ( ' eventHeader .mat ' ) ;  %Make  the  name  for  the  header 

.mat  file 


.mat  file 


wvf name  =  sprintf ( ' eventWaveform.mat ' ) ; %Make  the  name  for  the  waveform 


current  %channel 
stead  of  %period 
later  code) 


save  eventWaveform  waveform;  %Save  the  waveform 

traceName  =  char (traceNames (m) ) ;  %Make  a  string  out  of  the 

name 

samplingRate  =  1 /period (m) ;  %Make  a  sampling  rate  in¬ 

information 

% (more  useful  in 

save  eventHeader  traceName  samplingRate  shotDate  time;  %Save  the  header 
clear  waveform;  %Delete  the  waveform  variable  to  save 


memory 
event  time 


cd  .  . 


%Go  back  to  the  root  directory  for  the 


end 

end 

cd(rootDir)  %Back  to  the  day's  root  directory 

save  shotTimes  shotTimeArray;  %Save  the  list  of  event  times 

save  channels  channelArray;  %Save  the  list  of  channels  for  the  events 

fclose (wvf ID) ;  %close  the  .WVF  file 

f close (hdrID) ;  %close  the  .HDR  file 

end 


end 


end 
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clear 

pack; 

return; 
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function  [endian,  dataFormat,  dataOffset,  traceTotalNumber ,  traceNames,  blockSize,  vReso- 
lution,  vOffset,  vDataType,  period,  date,  time]  =  parseHeader ( f ile) 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Loads  the  header  .HDR  file  for  Yokogawa  DL750  waveform  objects,  and 
%extracts  necessary  data  out  of  it. 

%Written  for  experiments  in  Edgewood,  MD  June  2004  and  in  Blossom  Point,  MD 
%August  2004. 

%Version  1.0  completed  24  JUN  2004  by: 

%Tim  Eggerding 

% (t-eggerding@cecer.army.mil) 

%Update  history: 

%13  JUL  2004  -  added  "thisTrace"  index  in  order  to  handle  cases  where  a 
%  group  is  not  full. 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Referring  functions:  organizeWaveforms .m 


%Required  functions:  none 


%global  variables:  none 


%Inputs:  file 


File  name  for  the  header 


%Outputs:  endian 

%  dataFormat 


ming 
to  do 
the 

eluded  for 


%  dataOffset 

tion  takes  %up. 

%  traceTotalNumber 

%  traceNames 

%  blockSize 

%  vResolution 

%  vOffset 

%  vDataType 

bytes) 

%  period 

%  date 

%  time 


Holds  whether  big-endian  or  little-endian 
format 

"Trace"  or  "Block".  "Trace"  is  assumed 

througout  the  rest  of  the  code,  because  "Block"  program- 
does  not  make  sense  with  the  gathered  data.  "Block"  has 
with  saving  data  over  a  long  period  of  time,  which  is  not 
concern  of  this  study.  However,  this  variable  is  in¬ 
completeness  . 

How  many  bytes  within  the  .WVF  file  the  header  informa- 

How  many  channel's  worth  of  data  is  saved  in  this  header. 
Array  of  the  channel  names 

How  many  data  points  each  wave  consists  of 
The  resolution  of  the  screen  in  volts. 

Displacement  of  the  blast  wave  from  0  volts 

Encoding  of  the  how  the  data  points  should 

be  interpreted  (integer/flop,  signed/unsigned,  number  of 

The  period  of  sampling 
Date  of  event 
Time  of  event 


All  data  that  contains  textual  information  is  recorded  as  a  cell  array 
of  strings,  while  numbers  are  returned  as  a  numeric  array. 


%  See  the  DL750  operating  manual  page  APP-9  for  more  information  on  the 
%  structure  of  the  .HDR  file. 

%  Also  see  the  Yokogawa  Technical  Information  article  TI  7000-21E 
%  entitled  "Understanding  the  structure  of  Binary  data  (xxxxxxxx . WVF) 

%  file  created  by  the  DL,  AR  series."  A  copy  of  this  article  should  be 

%  with  the  DL750  manual  at  CER1,  but  can  also  be  found  at: 

%  http : //www. yokogawa . com/tm/pdf /tutorial/tm-technical_dlar . pdf 

headerText  =  textread ( f ile,  '%s');  %Read  in  the  header  file  as  a  cell  array  of 

strings 

endian  =  char (headerText ( 11 )) ;  %Store  the  endian  information  as  a  string 

dataFormat  =  char (headerText ( 13) ) ;  %Store  the  data  format  ast  a  string 

groupNumber  =  str2num (char (headerText ( 15) )) ;  %Store  the  number  of  groups 

traceTotalNumber  =  str2num (char (headerText ( 17 ))) ;  %Store  the  total  number  of  channels 

dataOffset  =  str2num (char (headerText (19) )) ;  %Store  the  number  of  bytes  of  header 

%information 
%index  variables 
%  i :  current  group 

%  j :  current  position  in  headerText  string  array 
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%  k:  current  trace  number  within  group 

j  =  22;  %j  points  to  the  point  in  the  headerText  array  that  holds  the  current  group 

name 

%For  loop  to  iterate  through  groups 

thisTrace  =  1; 

for  i  =  1 : groupNumber 

traceNumber  =  str2num (char (headerText (j ))) ;  %The  number  of  traces  (channels)  in 

the  group 

blockNumber  =  str2num (char (headerText (j+2) )) ;  %Number  of  blocks  per  trace.  For  the 
MD  %experiments, 

%this  should  always  be  1.  See  the  manual  for  the  dif¬ 
ference  b/w 

%Block  and  Trace  settings. 

%For  loop  to  iterate  through  traces  within  the  group 
for  k  =  1 : traceNumber 

traceNames (thisTrace)  =  headerText (j +3+k) ;  %Name  of  the  current  trace 
blockSize (thisTrace)  =  str2num (char (headerText ( j+3+traceNumber+l+k) ) ) ; 

%Number  of  data  points  saved  in  this  trace's  waveform 
vResolution (thisTrace)  =  str2num (char (headerText ( j+3+2* ( traceNumber+1 ) +k) ) ) ; 
%Resolution  of  the  data  in  volts 

vOf f set (thisTrace)  =  str2num (char (headerText ( j  +3+3* (traceNumber+1 ) +k) ) ) ; 

%DC  offset  of  the  data 

vDataType (thisTrace)  =  headerText (j+3+4* (traceNumber+1) +k) ; 

%Code  to  describe  the  size  and  format  of  data  points 
period (thisTrace)  =  str2num (char (headerText (j+3+11* (traceNumber+1) +k) ) )  ; 

%Sampling  period 

date (thisTrace)  =  headerText (j+3+14* (traceNumber+1) +k) ; 

%Date  of  event 

time (thisTrace)  =  headerText (j+3+15* (traceNumber+1) +1) ; 

%Time  of  event 
thisTrace  =  thisTrace+1; 

end 

j=j+3+16* (traceNumber+1) +2;  %Increment  j  to  point  at  the  next  group 

end 
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function  writeSensitivities (calLevel,  calFreq, directory) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Loop  through  a  directory  of  calibration  signals  and  write  the  sensitivies 
%given  by  f indSensitivity 

%Written  for  experiments  in  Edgewood,  MD  June  2004  and  in  Blossom  Point,  MD 
%August  2004. 

%Version  1.0  completed  24  JUN  2004  by: 

%Tim  Eggerding 

% (t-eggerding@cecer.army.mil) 

%Update  history:  none. 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 


%Ref erring  functions:  yokogawaDataGUI .m 


%Required  functions: 


f indSensitivity .m 
parseHeader .m 


%global  variables:  none 


%Inputs:  calLevel 

%  calFreq 

%  directory 


SPL  of  calibrator  in  dB 
Frequency  of  calibrator  in  Hz 
Directory  containing  calibration  files 


%Outputs:  sensDat.mat 


Calibration  in  V/Pa 


cd (directory) ; 

D  =  dir; 

numFiles  =  length (D); 
for  i  =  1: numFiles 
wvfName  =  D(i) .name; 


%Creates  a  vector  listing  all  the  files  in  the  directory 

%numFiles  =  the  number  of  files  in  the  directory 

%Loop  through  the  files 

%Take  the  name  of  the  current  file 


if (~isempty (strf ind (wvfName, ' WVF ' ) ) ) 
hdrName  =  [wvfName ( 1 : 4 )  ' . HDR']; 

hdrID  =  f open (hdrName,  '  r'); 
if  hdrID  ==  -1 


%If  this  is  a  .WVF  file  (thus,  a  calibration) 
%Create  the  name  of  the  appropriate  header  file 
%Open  the  header  file 


%Warn  if  the  header  is  not  found 

warning ([' Header  file  '  hdrName  '  for  waveform  file'  wvfName  '  does  not  ex¬ 
ist.']);  %display  error  message  if  no  header  file, 
else 


%Get  the  necessary  data  from  the  header  file 

%In  the  case  of  a  calibration,  variables  like  traceNames,  which 
%are  usually  arrays,  will  only  have  one  value,  since  only  one  cal 
%is  saved  at  a  time. 

[endian,  dataFormat,  dataOffset,  traceTotalNumber ,  traceNames,  blockSize, 
vResolution,  vOffset,  vDataType,  period,  date,  time]  =  parseHeader (hdrName) ; 

traceNames  =  char (traceNames)  %Cast  traceNames  as 


string 

string 

sary 


vDataType  =  char (vDataType) ; 

if  strcmp (endian,  'Big') 

wvfID  =  f open (wvfName, ' r ' ,  'ieee-be'); 

else 

wvfID  =  fopen  (wvfName,  'r',  'ieee-le'); 


%Cast  vDataType  as 

%Big-endian  if 

%Little-endian 


a 

a 


neces- 


if 


necessary 

end 

headerData  =  f read (wvfID, dataOff set/2 ,  ['*intl6']);  %Read  the  header 

data 


clear  headerData 
%and  throw  it  out 


for  m=l : size (traceNames, 1 ) 

if (~ (exist (traceNames (m, :) ,  'dir')  ==  7)) 
for  this  %channel  has  not  occured 


channel 


mkdir (traceNames (m, : ) ) ; 
cd (traceNames (m, : ) ) ; 


%If  a  calibration 

%Make  a  directory  for  this 
%Enter  that  di¬ 


rectory 


waveform  =  fread (wvfID,  blockSize (m) ,  ['*intl6']); 
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sensDat  =  findSensitivity (double (waveform) *vResolution (m) , 
1/period (m),  calFreq,  calLevel) ; 

filename  =  [' sensDat .mat '] ; 
save (filename,  'sensDat'); 
cd  .  . 


end 

end 

end 

end 

end 

fclose ( ' all '  ) 

cd  .  . 
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function  [sensitivity] =f indSensitivity  (data, samplingRate,  calFreq,  calLevel) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Takes  the  FFT  of  the  data  and  eliminates  the  high  and  low  frequency  data 
%of  the  calibration  signal.  The  inverse  FFT  is  then  taken  and  the 
%sensitivity  is  found  in  V/Pa  by  Vrms/Prms. 

%Written  for  experiments  in  Edgewood,  MD  June  2004  and  in  Blossom  Point,  MD 
%August  2004. 

%Adopted  from  code  written  by  Ryan  Lee  for  the  Texarkana  experiments. 

%Version  1.0  completed  24  JUN  2004  by: 

%Tim  Eggerding 

% (t-eggerding@cecer.army.mil) 

%Update  history:  none. 

^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 


%Referring  functions:  writeSensitivity .m 


%Required  functions:  none 

%global  variables:  none 


%Inputs : 


data 

samplingRate 

calFreq 

calLevel 


Calibration  signal  waveform  given  in  Volts 
Rate  of  sampling  from  the  scope 
Frequency  of  the  calibrator  used  in  Hz 
SPL  of  the  calibrator  used  in  dB 


%Outputs:  sensitivity 


Calibration  in  V/Pa 


pO  =  20*10A-6; 

lengthData  =  length (data) ; 

lengthFFT  =  2 A (nextpow2 ( length (data) ) -1 ) ; 

be  %truncated  (rather 


freqRes  =  samplingRate/lengthFFT; 
domain 


%Reference  pressure  in  Pascals 

%This  will  be  the  length  of  the  FFT 
%and  the  "-1"  insures  that  our  data  will 

%than  zero  padded)  in  the  FFT  algorithm 

%Resolution  of  frequencies  in  the  frequency 


H  =  fft  (data, lengthFFT) ;  %  This  is  the  FFT  in  "H" 

%  IMPORTANT  NOTE:  H(l)  is  really  the  element  of  the 
%  FFT  that  is  at  a  frequency  of  0 .  Matlab  does  not 
%  index  any  matricies  beginning  with  0,  it  always  be 


gins 

%  with  an 
%  Matlab. 

cutoff_l  =  round (calFreq*2A (-1/12 ) /freqRes) ; 
eliminate  %noise 

cutoff_2  =  round (calFreq*2A ( 1/12 ) /freqRes) ; 

cutoff_3  =  length (H)  -  cutoff_2; 

cutoff_l,  %since 

cutoff_4  =  length (H)  -  cutoff_l; 

are  %actually 


ndex  of  1.  Therefore  H(0)  is  improper  in 

%Calculate  cutoff  frequencies  to 

%calFreq  +/-  a  semitone 

%The  negative  cutoff  frequency  of 

%the  values  of  a  DFT  greater  than  pi 

%negative  frequencies 


H ( 1 : cutof f_l )  =  0; 


%Zero  out  all  unwanted  low  and  high  frequencies  from 
%zero  to  one  semitone  below  our  calibration  frequency 


one 

H (cutof f  2 : cutof f  3)  =  0; 


%Zero  from  one  semitone  above  our  calibration  frequency  to 
%semitone 

%below  the  negative  calibration  frequency 


%Zero  From  one  semitone  above  the  negative  calibration  fre¬ 
quency  to 

H (cutof f_4 : length (H) )  =  0;  %the  end  of  the  FFT. 

data  =  if ft (H, length (H) ) ;  %New  time-based  signal.  This  should  be  the  original 
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zero . 

Vrms  =  sqrt (sum (abs (data . A2 
the  mean 

Prms  =  10 A (calLevel/20 ) *p0 ; 
SPL 

sensitivity  =  -Vrms/Prms; 
return 


%cal  signal  minus  any  unwanted  noise.  Its  mean  should  be 

)) /length (data) ) ;  %RMS  voltage  found  by  the  square  root  of 

%of  the  square  of  the  data. 

%RMS  pressure  found  by  the  inverse  of  the 

%equation 

%SPL(dB)  =  20*logl0 (Prms/POrms) 
%Microphone  sensitivity  in  V/Pa. 
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function  plotonl (rootdir) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  directory  for  one  shot 
%Output :  none 

%Description :  Loads  each  channel's  waveform  for  a  single  shot  and  plots  the  calibrated 
%signals  all  on  one  graph  in  Pa  over  seconds. 

%Requires :  eventWaveform.mat , eventHeader .mat ,  sensDat .mat , MicLocations .mat 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 


cd  (rootdir) ; 

D  =  dir (rootdir)  ; 
numFiles  =  length (D); 
loc=cell (numFiles-2 ,  1)  ; 
count=0 ; 


%Takes  as  input  the  directory  of  one  shot 

%Finds  how  many  channels  in  the  directory 
%Preallocate  space  for  legend  array 
%Initialize  count  variable 


for  i  =  1: (numFiles) 

channel= [ ' CH '  num2str(i)]; 


%Scrolling  through  all  channels 


if  exist (channel, ' dir ') ~=7  %Checks  if  there's  data  for  the  channel 
continue 

end 


count=count+l ; 

cd  (channel)  %0pens  folder  for  each  channel 

load  eventWaveform;  %and  loads  waveform  in  Volts/Sample 

load  eventHeader;  %as  well  as  samplingRate,  shotDate,  and  time 

cd  (['../.. /PRECAL/ '  channel] ) 

load  sensDat;  %Gets  sensitivity  data  from  calibration 


%loads  mic  locations  for  ei- 
%forest  scope  or  the  trailer 
%display  in  legend 

elseif  length (findstr (rootdir,  ' Trailer ')) ~=  0 
loc { count, 1}  =  TrailerMicLocations { i, 1 } ; 

end 


cd  (['../../..'])  ; 
ther  the 

load  MicLocations .mat; 

for 

if  length (findstr (rootdir,  'Forest' ))~=  0 
loc {count,!}  =  ForestMicLocations { i, 1 } ; 


timeaxis=[l: (length (waveform) ) ] '/samplingRate; 
time  (s) 

calwaveform  =  waveform/ sensDat; 
if  (mod (i, 2) ==0) 

plot  (timeaxis, calwaveform,  'r') 

else 

plot  (timeaxis, calwaveform,  'b') 

end 


%Sets  up  the  X-axis  to  show 

%Converts  waveform  to  Pascals 

%and  plots  calibrated  waveform 
%alternately  in  red  and  blue 


hold  on  %holds  plot  so  others  will  plot  on  top  of  it 

cd  (rootdir) 

end 


hold  off 

xlabel('Time  (s)'); 
ylabel (' Pressure  (Pa)'); 

title ([' Sample  from  '  shotDate  '  at  '  time  '.']); 
legend (loc ( 1 : count) , ' Location ' , ' NorthEastOutside ' ) ; 
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function  plotsingle (rootdir,  chanNum) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  directory  for  one  shot,  channel  number 
%Output :  none 

%Description :  This  function  will  take  as  input  the  directory  for  a  single  shot  and  the 
%channel  whose  plot  is  desired.  It  will  then  plot  the  calibrated  waveform  in 
%Pascals.  It  will  also  indicate  the  peak  level  of  the  waveform. 

%Requires :  eventWaveform.mat, eventHeader .mat, sensDat .mat, MicLocations .mat 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 


channel  =  ['CH'  num2str (chanNum) ] ; 

cd  (channel) 

load  eventHeader .mat; 

time 

load  eventWaveform.mat; 

cd  (['../.. /PRECAL/ '  channel] ) ; 
load  sensDat. mat; 

timeaxis=[l: (length (waveform) ) ] ' /samplingRate; 
(sec) 

calwaveform  =  waveform/ sensDat; 
peakLevel  =  max (abs (calwaveform) ) ; 
neg) 

cd (['../../..'])  ; 

the 

load  MicLocations .mat ; 

if  length (findstr (rootdir,  'Forest' ))~=  0 
loc  =  ForestMicLocations { chanNum, 1 } ; 

elseif  length (findstr (rootdir,  ' Trailer ')) ~=  0 
loc  =  TrailerMicLocations { chanNum, 1 } ; 

end 

plot (timeaxis, calwaveform) 
xlabel('Time  ( s ) ' ) ; 
ylabel (' Pressure  (Pa)'); 

title ([' Sample  from  '  shotDate  '  at  '  time  '. 
Pa.']); 

legend (loc, ' Location ' , ' NorthEastOutside ' ) ; 


%Load  samplingRate,  shotDate, 
%and  the  uncalibrated  waveform 

%Load  sensitivity  data 

%Sets  up  the  X-axis  to  show  time 

%Converts  waveform  to  Pascals 
%Calculates  peak  level  (pos  or 

%loads  mic  locations  for  either 

%forest  scope  or  the  trailer  for 
%display  in  legend 


%Plots  calibrated  waveform 

Peak  Level  =  '  num2str (peakLevel)  ' 


cd  (rootdir) 
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function  plotPeaks  (srcin,  srcloc,  config) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  (all  integers)  source  type,  source  location,  array  config.  Based  on  the  fol- 


lowing 

table : 

% 

SOURCE 

LOCATION 

CONFIG 

% 

1  =  C-4 

1  =  BP1 

1  =  Parallel 

% 

2  =  PC 

2  =  BP2 

2  =  Perpendicular 

% 

3  =  BP3 

3  =  Vertical 

%0uput : 

:  none 

4  =  BP4 

4  =  All 

%Description :  Creates  one  subplot  that  shows  the  peak  level  for  every  mic 

%in  dB  over  distance  from  source.  A  second  subplot  shows  a  horizontal  cross 
%section  of  the  field  (with  the  forest  edges  indicated)  and  shows  the  peak 
%levels  at  each  mic  location. 

%Requires:  peaksTable .mat 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

load  ('N:\Research\Acoustics\Camp  Ripley  Matlab  Format\Pk  vs  Dist  Analy- 
sis\peaksTable .mat ' ) ; 
elf  ; 


Al=zeros (length (peaksTable) 
loc  where 

,2) 

;  Alcnt=0; 

%Initializes  a  matrix 

for  each  mic 

A2=zeros (length (peaksTable) 

,2) 

;  A2cnt=0; 

%col  1  =  dist  from  blast  point  and 

col  2  = 

A3=zeros (length (peaksTable) 

,2) 

;  A3cnt=0; 

%SPL  in  dB 

A4=zeros (length (peaksTable) 

,2) 

;  A4cnt=0; 

Ll=zeros (length (peaksTable) 
plotting 

,2) 

;  Llcnt=0; 

%The  *cnt  variable  is 

used  later  in 

L2=zeros (length (peaksTable) 

,2) 

;  L2cnt=0; 

%to  keep  from  running 

into  errors 

L3=zeros (length (peaksTable) 

,2) 

;  L3cnt=0; 

Sl=zeros (length (peaksTable) 

,2) 

;  Slcnt=0; 

S2=zeros (length (peaksTable) 

,2) 

;  S2cnt=0; 

S3=zeros (length (peaksTable) 

,2) 

;  S3cnt=0; 

S4=zeros (length (peaksTable) 

,2) 

;  S4cnt=0 ; 

S5=zeros (length (peaksTable) 

,2) 

;  S5cnt=0; 

for  i=l : length (peaksTable) 

%Goes  through  table  and  finds  all  rows  that  match  the  desired  source  type, 

%source  location,  and  either  a  particular  config  or  all  configs.  The  dist 
%from  the  source  and  the  SPL  are  stored  in  the  corresponding  matrix  for 
%that  mic. 

%peaksTable  is  set  up  with  columns: 

%  l=Source  type  2=Source  Loc  3=Array  Config  4=Mic  Loc  5=dist  from  source 
%  6=SPL  in  Pa  7=SPL  in  dB 

if  (peaksTable (i, 1) ==srcin  &&  peaksTable (i, 2) ==srcloc  &&  (con- 
fig==peaksTable (i, 3) I |config==4)) 

if  (ltpeaksTable (i, 4)  &&  peaksTable (i, 4) <=4) 

Alcnt=Alcnt+l ; 

A1 (Alcnt , 1 ) =peaksTable ( i , 5 ) ;  A1 (Alcnt , 2 ) =peaksTable ( i , 7 ) ; 
elseif  (4<peaksTable (i, 4)  &&  peaksTable (i, 4) <=8) 

A2cnt=A2cnt+l ; 

A2 (A2cnt , 1 ) =peaksTable ( i , 5 ) ;  A2 (A2cnt , 2 ) =peaksTable ( i , 7 ) ; 
elseif  (8<peaksTable (i, 4)  &&  peaksTable (i, 4) <=12) 

A3cnt=A3cnt+l ; 

A3 (A3cnt, 1) =peaksTable (i, 5) ;  A3 (A3cnt, 2) =peaksTable (i, 7) ; 
elseif  (12<peaksTable (i, 4)  &&  peaksTable (i, 4) <=16) 

A4cnt=A4cnt+l ; 

A4 (A4cnt , 1 ) =peaksTable ( i , 5 ) ;  A4 (A4cnt , 2 ) =peaksTable ( i , 7 ) ; 
elseif  (peaksTable (i, 4) ==17) 

Llcnt=Llcnt+l ; 

LI (Llcnt , 1 ) =peaksTable ( i , 5 ) ;  LI (Llcnt , 2 ) =peaksTable ( i , 7 ) ; 
elseif  (peaksTable (i, 4) ==18) 

L2cnt=L2cnt+l ; 
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L2 (L2cnt , 1 ) =peaksTable ( i , 5 ) ;  L2 (L2cnt , 2 ) =peaksTable ( i , 7 ) ; 
elseif  (peaksTable (i, 4) ==19) 

L3cnt=L3cnt+l ; 

L3 (L3cnt, 1) =peaksTable (i, 5) ;  L3 (L3cnt, 2) =peaksTable (i, 7) ; 
elseif  (peaksTable (i, 4) ==20) 

Slcnt=Slcnt+l ; 

SI (Slcnt, 1) =peaksTable (i, 5) ;  SI (Slcnt, 2) =peaksTable (i, 7) ; 
elseif  (peaksTable (i, 4) ==21) 

S2cnt=S2cnt+l ; 

S2 (S2cnt, 1) =peaksTable (i, 5) ;  S2 (S2cnt, 2) =peaksTable (i, 7) ; 
elseif  (peaksTable (i, 4) ==22) 

S3cnt=S3cnt+l ; 

S3 (S3cnt, 1) =peaksTable (i, 5) ;  S3 (S3cnt, 2) =peaksTable  (i, 7) ; 
elseif  (peaksTable (i, 4) ==23) 

S4cnt=S4cnt+l ; 

S4 (S4cnt, 1) =peaksTable (i, 5) ;  S4 (S4cnt, 2) =peaksTable (i, 7) ; 
elseif  (peaksTable (i, 4) ==24) 

S5cnt=S5cnt+l ; 

S5 (S5cnt, 1) =peaksTable (i, 5) ;  S5 (S5cnt, 2) =peaksTable (i, 7) ; 

end 

end 

end 

%First  subplot  shows  level  vs  actual  distance  from  the  blast  point 
subplot (2,1,1) 

plot (A1 ( 1 : Alcnt , 1 ) , A1 ( 1 : Alcnt , 2 ) , ' bo ' ) 
hold  on 

plot (A2 ( 1 : A2cnt , 1 ) , A2 ( 1 : A2cnt , 2 ) , ' bs ' ) 
plot (A3 ( 1 : A3cnt, 1 ) , A3 ( 1 : A3cnt, 2 ) , ' bd ' ) 
plot (A4 ( 1 : A4cnt , 1 ) , A4 ( 1 : A4cnt , 2 ) , ' bp ' ) 
plot (LI ( 1 : Llcnt , 1 ) , LI (1 : Llcnt , 2 ) , ' mo ' ) 
plot (L2 ( 1 : L2cnt , 1 ) , L2 ( 1 : L2cnt , 2 ) , ' ms ' ) 
plot (L3 ( 1 : L3cnt, 1 ) , L3 (1 : L3cnt, 2 ) , 'md ' ) 
plot (SI (1 : Slcnt, 1) , SI (1 : Slcnt, 2) , ' go ' ) 
plot (S2 (1 :S2cnt, 1) , S2 (1 :S2cnt,2) , ' gs ' ) 
plot (S3 (1 :S3cnt, 1) , S3 (1 :S3cnt,2) , 'gd' ) 
plot (S4 (1 :S4cnt, 1) , S4 (1 : S4cnt,2) , ' gp ' ) 
plot (S5 ( 1 : S5cnt, 1 ) , S5 (1 : S5cnt, 2 ) , ' g+ ' ) 

grid  on 

legend ( ' A1 1 ,  ' A2 ' ,  ' A3  1 ,  ' A4 ' , 1  LI '  ,  1 L2 ' , 1 L3 ' ,  1  SI 1 ,  ' S2  '  ,  '  S3  '  ,  '  S4  1  ,  '  S5  '  ,  .  .  . 

' Location '  ,  ' NorthEastOutside ' ) 
xlabel (' Distance  from  Source  (meters)') 
y label ( 1 SPL  ( dB ) ' ) 
if  srcin==l 

source= ' C-4 ' ; 
elseif  srcin==2 
source= 'PC'; 

end 

if  config==l 

arrayConf ig= ' Parallel ' ; 
elseif  config==2 

arrayConf ig= ' Perpendicular ' ; 
elseif  config==3 

arrayConf ig= ' Vertical ' ; 
elseif  config==4 

arrayConf ig= ' All ' ; 

end 

title ([ source  '  -  BP'  num2str (srcloc)  '  :  '  arrayConfig  '  Config']) 

hold  off 

%Second  subplot  shows  a  horizontal  view  where  the  mic  locations  are 
%stationary  and  the  blast  point  moves  along  the  horizontal  axis 

BPlxy=[0  0;0  200];  %Define  x  and  y  coordinates  for  the  blast  points  and  the  forest  edges 
BP2xy= [ 0  0 ; 0  200]; 

BP3xy= [175  0; 175  200]  ; 

BP4xy= [313  0;  313  200]; 
edgel= [50  0;  50  200]; 
edge2= [263  0;  263  200]; 

%Define  x  and  y  coordinates  for  mic  Iocs  along  main  prop  line 
Alxy= [0  25] ; 
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A2xy= [ 0 

50]  ; 

A3xy= [0 

75]  ; 

A4xy= [0 

100]  ; 

Llxy= [ 0 

150]  ; 

L2xy= [ 0 

2  63]  ; 

L3xy= [0 

288]  ; 

subplot 

(2,1,2 

plot (edgel {:,!), edgel  ( : , 2 ) , 
hold  on 

k- 

-')  ; 

plot (edge2 {:,!), edge2 ( : , 2 )  , 
if  srcloc==l 

k- 

-1)  ; 

plot (BPlxy ( : , 1 ) , BPlxy ( : 
elseif  srcloc==2 

2) 

,  '  r- 

plot (BP2xy {:,!), BP2xy ( : 
elseif  srcloc==3 

2) 

,  '  r- 

plot (BP3xy {:,!), BP3xy ( : 
elseif  srcloc==4 

2) 

,  '  r- 

plot (BP4xy ( : , 1 ) , BP4xy ( : 

2) 

,  '  r- 

end 

legend  ( ' Edge ' , ' Edge '  ,  ' BP ' ,  . . . 

' Location ' , ' NorthEastOutside ' ) ; 
stem (Alxy ( 1 , 2 ) *ones (Alcnt, 1 ) ,  A1 (1 : Alcnt, 2 )  , ' ob ' ) ; 
stem (A2xy (1,2) *ones (A2cnt, 1 ) ,A2  (1 : A2cnt, 2) , ' sb ' ) ; 
stem (A3xy (1,2) *ones (A3cnt, 1) , A3 (1 : A3cnt, 2) , ' db ' ) ; 
stem (A4xy ( 1 , 2 ) *ones (A4cnt, 1 ) , A4 (1 : A4cnt, 2 ) , ' pb ' ) ; 
stem (Llxy ( 1 , 2 ) *ones (Llcnt , 1 ) , LI ( 1 : Llcnt , 2 ) , ' om ' ) ; 
stem (L2xy ( 1 , 2 ) *ones (L2cnt, 1 )  ,  L2 (1 : L2cnt,  2 )  ,  ' sm' ) ; 
stem (L3xy (1,2) *ones (L3cnt , 1 ) , L3 (1 : L3cnt , 2 ) , ' dm ' ) ; 

hold  off 

xlabel (' Horizontal  Distance  (meters) '); ylabel (' SPL  (dB) ' ) ; 
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function  peaksOut=getPeaks (rootdir) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  directory  for  one  day's  sample  set 

%Output :  matrix  indicating  shot  date,  shot  time,  channel  number,  and  peak  levels 

%for  each  waveform  from  that  day 

%Description :  For  use  in  the  Peak  level  vs  Distance  Analysis: 

%Takes  as  input  the  directory  for  one  day's  samples,  and  returns  a  matrix 
%indicating  the  peak  levels  for  each  waveform  from  that  day. 

%Requires :  eventWaveform.mat, eventHeader .mat, sensDat .mat 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

cd (rootdir) 

D=dir (rootdir) ;  numf iles=length (D) ; 
count=0;  peaks=cell ( (numf iles-2 ) *  1 6 ,  4); 

for  j  =  lznumfiles 

if  ( D  ( j  ,  1 )  .  isdir==l  &&  strcmp  (D  ( j  ,  1 )  .  name,  '.')==()  &&  strcmp  ( D  ( j  ,  1 )  .  name,  '  PRECAL  1  )  ==0 
&&  strcmp (D ( j , 1 ) . name, ' POSTCAL ' ) ==0) 

cd (D ( j , 1 ) . name) 

for  i  =  1:16  %Scrolling  through  all  channels 

channel= [ ' CH '  num2str(i)]; 

if  exist ( channel ,' dir ') ~=7  %Checks  if  there's  data  for  the  channel 
continue 

end 

count=count+l ; 

cd  (channel)  %0pens  folder  for  each  channel 

load  eventWavef orm;  %and  loads  waveform  in  Volts/Sample 

load  eventHeader;  %as  well  as  samplingRate,  shotDate,  and  time 

cd  (['../../PRECAL/'  channel]) 

load  sensDat;  %Gets  sensitivity  data  from  calibration 

calwaveform  =  waveform/sensDat;  %Converts  waveform  to 

Pascals 

peaks {count,  1}  =  shotDate; 
peaks {count,  2}  =  time; 
peaks{count,  3}  =  traceName; 
peaks{count,  4}  =  max (abs (calwaveform) ) ; 

cd (['../ _ /'  D ( j , 1 ) . name ] ) 

end 

end 

cd  (rootdir) 

end 

peaksOut=peaks ( 1 : count  ,  1:4); 
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function  WaveandSpectra  (rootdir,  chanNum) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  directory  for  one  shot,  channel  number  of  desired  mic 

%Output :  none 

%Description :  This  function  will  take  as  input  the  directory  for  a  single  shot  and  the 

%channel  whose  plot  is  desired.  It  will  then  plot  the  calibrated  waveform  in 
%Pascals  and  the  1/3  octave  band  spectra.  Default  beginning  and  end  times  for  the  FFT  are 
0.7s 

%and  3.7s  respectively.  It  will  also  indicate  the  peak  level  of  the  waveform. 

%Requires:  eventHeader .mat ,  eventWaveform.mat ,  sensDat,  MicLocations .mat 

%Calls:  SELtest2.m,  semilogxBar .m 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

channel  =  [ ' CH '  num2str (chanNum) ] ; 
cd  (channel) 

load  eventHeader .mat;  %Load  samplingRate,  shotDate, 

time 

load  eventWaveform.mat;  %and  the  uncalibrated  waveform 


cd  (['../.. /PRECAL/ '  channel]); 
load  sensDat .mat; 

timeaxis=[l: ( length (waveform) ) ] '/samplingRate; 
(sec) 

calwaveform  =  waveform/ sensDat; 
peakLevel  =  max (abs (calwaveform) ) ; 
neg) 

cd (['../../..'])  ; 

the 

load  MicLocations .mat; 

if  length (findstr (rootdir,  'Forest' ))~=  0 
loc  =  ForestMicLocations { chanNum, 1 } ; 

end 

if  length (findstr (rootdir,  ' Trailer ')) ~=  0 
loc  =  TrailerMicLocations { chanNum, 1 } ; 

end 

subplot  (2,1,1),  plot ( timeaxis , calwaveform) 
xlabel('Time  ( s ) ' ) ; 
ylabel ( ' Pressure  (Pa)  '  )  ; 

title ([' Sample  from  '  shotDate  '  at  '  time  '. 
Pa .  '  ]  )  ; 

legend (loc,  ' Location ' , ' Nor thEastOut side ' ) ; 

tl=0.7;  %define  start  and  stop  times  i 

t2=3 . 7 ; 

[tobSEL,  Fcenter,  totalSEL]  =  SELtest2 (rootdii 


%Load  sensitivity  data 

%Sets  up  the  X-axis  to  show  time 

%Converts  waveform  to  Pascals 
%Calculates  peak  level  (pos  or 

%loads  mic  locations  for  either 

%forest  scope  or  the  trailer  for 
%display  in  legend 


%and  plots  calibrated  waveform 
Peak  Level  =  '  num2str (peakLevel)  ' 

or  FFT 

,  chanNum,  tl,  t2); 


subplot ( 2 , 1 , 2 ) ,  semilogxBar (Fcenter,  tobSEL);  %Generate  and  plot  the  sound  ex¬ 

posure  %level 
xlabel ( ' Frequency  (Hz) ' ) ; 
ylabel ( ' SEL  (dB)  ' )  ; 

axis ( [sqrt (10A (-2/5) ) ,  sqrt ( 1 0  A ( (length (Fcenter) +1) /5) ) ,  10*floor(. l*min (tobSEL) ) , 
10*ceil (max (tobSEL) * . 1 ) ] ) 


cd  (rootdir) 
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function  [tobSEL, Fcenter , totalSEL] =SELtest2 (rootdir,  chanNum, tl , t2 ) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  directory  for  one  shot,  desired  channel,  start  time,  stop  time 
%0utput:  array  of  the  1/3  octave  band  SEL  levels,  center  freqs  of  bands 
%Description :  Calculates  and  returns  the  1/3  octave  band  SEL  for  one  shot 
%waveform.  The  user  defines  the  beginning  and  ending  times  to  be  considered 
%for  this  calculation,  where  time  is  from  0  to  5  sec. 

%Requires:  eventWavef orm . mat ,  eventHeader .mat ,  sensDat.mat 
%Calls:  generateBands .m,  semilogxBar .m 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

cd  (rootdir) 

channel  =  ['CH*  num2str (chanNum) ] ; 

%If  there  is  no  data  for  that  channel,  return  zeros  for  all  outputs 
if  exist (channel, ' dir ') ==0 

tobSEL= zeros (1,47) ; Fcenter= zeros (1,47) ; total SEL=0 ; 
return 

end 

cd  (channel) 
load  eventHeader .mat; 
time 

load  eventWavef orm. mat; 

cd  (['../.. /PRECAL/ '  channel] ) ; 
load  sensDat.mat;  %Load  sensitivity  data 

cd  (rootdir) 

if  tl>t2 

temp=tl;  %Makes  sure  tl<t2 

tl=t2 ; 

t2=temp; 

end 

calwaveform  =  waveform/sensDat ;  %Converts  waveform  to  Pascals 

%signal  will  be  the  segment  of  the  wave  we  want  to  evaluate 
signal=calwaveform (round (tl*samplingRate) +1 : round (t2*samplingRate) , 1 ) ; 

pO  =  20e-6;  %Define  reference  pressure 


%Load  samplingRate,  shotDate, 
%and  the  uncalibrated  waveform 


S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Get  1/3  Oct  Band  SEL  for  signal 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

N  =  2 Anextpow2 (samplingRate* (t2-tl )) ;  %Find  the  size  of  the  FFT  necessary  for  t2-tl  sec¬ 
onds  of  %data 

freqRes  =  samplingRate/N;  %The  resolution  of  frequencies  in  the  frequency  domain 

lowRelFreq  =  floor (samplingRate/N) ;  %Lowest  reliable  frequency  in  the  FFT 

%Generate  bands  for  the  SEL  using  generateBands .m 

[Flo, Fcenter, Fhi]  =  generateBands (lowRelFreq, samplingRate) ; 

%Calculate  FFT 
H  =  fft  (signal, N); 

%Compute  Narrow  Band  Sound  Exposure 
exposure  =  abs (H) . A2/ (samplingRate*N) ; 
totalSEL  =  10*logl0 (2*sum (exposure) /pOA2 ) ; 

%Calculate  1/3-Octave  Band  Sound  Exposure 

indexLow  =  round (Flo ( 1 ) /freqRes) ;  %The  first  frequencies  added  up  will  be 

%starting  at  the  lowest  Flo 

if  round (Fhi ( 1 ) /freqRes)  >  1  %Must  take  care  of  the  first  band  calcula¬ 

tions 

indexHigh  =  round (Fhi ( 1 ) /freqRes) -1 ;  %Separately  here,  to  avoid  indexing  by  0 
tobShot ( 1+lowRelFreq)  =  2*sum (exposure (indexLow: indexHigh) ) ; 
indexLow  =  indexHigh  +  1; 

end; 
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for  k=2 : length (Flo)  %  For  k  =  2:#  of  Bands,  for  the  rest  of  the  bands 

indexHigh  =  round (Fhi (k) /freqRes) -1 ; 
if  indexHigh  >=  indexLow 

if  round (Fhi (k) /freqRes)  >  round (Fhi (k-1 ) /freqRes) 

tobShot (k+lowRelFreq)  =  2*sum (exposure (indexLow: indexHigh) ) ;  %T.O.B.  is  one- 

Third 

%Octave  Band 

sound  exposure 

indexLow  =  indexHigh  +  1; 

end; 

end; 

end; 


%Compute  1 /3-Octave  Band  SEL 

tobSEL  =  10*logl0 (tobShot/pO A2 ) ;  %l/3  octave  band 

%sound  exposure  level  =  10  log  (E  i/P  oA2) 


ERDC  SR-07-3 


52 


function  [Fcenter, totalSEL, tobSEL] =SELavg2 (srcType, srcLoc, scope, channel) 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  srcType=  'C-4'  or  'PC' 

%  srcLoc=  1,2, 3, 4 

%  scope=  'Trailer'  or  'Forest'  (must  be  entered  exactly  like  these) 

%  channel=  1  thru  16 

%Output :  Fcenter=  47  element  vector  with  center  freqs  for  the  1/3  oct  bands 

%  totalSEL=  total  SEL  for  each  shot 

%  tobSEL=  each  row  represents  one  shot,  where  the  cols  are  the  SEL 

%  for  each  band  corresponding  to  the  same  col  in  Fcenter. 

%Description :  Loads  DataLog.mat  and  filters  to  gather  all  shots  matching 

%  the  inputs.  Then  runs  each  shot  through  SELtest2.m  to  get 

%  tob  SEL  values.  Levels  for  each  shot  are  compiled  into  tobSEL. 

%  Default  values  for  FFT  start  and  stop  times  are  0.7sec  & 

%  3.7sec.  Average  SEL  level  across  each  shot  are  also  computed 

%  and  can  be  returned  or  plotted. 

%Requires:  DataLog.mat 

%Calls:  SELtest2.m,  semilogxBar .m 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

cd  ('N:\Research\Acoustics\Camp  Ripley  Matlab  Format'); 
load  DataLog.mat; 
count=0 ; 
x=  [  ]  ; 

%DataLog  contains  information  for  every  shot.  Its  columns  are: 

%  1) Scope  2)Yok  file  number  3) Date  4) Time 

%  5) Source  type  6) Source  location  7) Array  config 

%  8) file  location (within  N:\Research\Acoustics\Camp  Ripley  Matlab  Format) 

%Shots  in  DataLog  that  match  the  desired  inputs  are  collected  into  array  x 
%which  will  have  length  of  count 
for  i=l : length (DataLog) 

if  strcmp (DataLog{ i, 6 } , [ ' BP '  num2str (srcLoc) ]) ==1  ... 

&&  strcmp (DataLog{ i, 5 } ,  srcType)==l  ... 

&&  strcmp (DataLog! i, 1 } ,  scope) ==1 
count=count+l ; 
for  z=l : 8 

x { count , z } =DataLog { i ,  z }  ; 

end 

end 

end 

if  isempty (x) ==1 

disp ( ' There  are  no  waveforms  that  match  the  input  values.') 

tobSEL=zeros (1,47) ; Fcenter=zeros (1,47) ; totalSEL=0 ; 

return 

end 

tl=0.7;  %Define  start  and  stop  times  for  the  FFT  calculation 

t2=3 . 7 ; 

%Fcenter  contains  the  center  frequencies  for  the  1/3  octave  bands 

%For  each  shot  in  x,  there  will  be  a  row  in  tobSEL  that  lists  the  amplitude 

%for  each  corresponding  frequency  band. 

%The  total  SEL  for  each  shot  will  be  collected  in  array  totalSEL 
for  k=l: count 

[tobSEL (k, : ) , Fcenter, totalSEL (k) ] =SELtest2 ( [ ' N : \Research\Acoustics\Camp  Ripley  Matlab 
Format'  x{k, 8}],  ... 

channel,  tl,t2); 

end 

%To  find  average  SEL  across  all  the  shots,  all  the  columns  in  tobSEL  are  summed 
%and  divided  by  the  number  of  shots 
for  a=l : length (tobSEL ' ) 

tobSELavg (a) =sum (tobSEL ( : , a) ) /count; 

end 

totalSELavg=sum (totalSEL) /length (totalSEL) ; 
semilogxBar (Fcenter ,  tobSELavg) 
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function  [results] =collectSEL2 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%Input:  None  (user  defines  srctype,  srcloc,  and  scope  within  m-file: 

%  line  23-26) 

%Output:  Matrix  containing  source  type,  source  location,  mic  location, 

%  total  SEL  (in  dB) ,  and  1/3  octave  band  SEL  levels  for  the  center  frequencies 

%  given  in  row  1  cols  5+ 

%Description :  For  the  user  defined  values,  the  function  runs  SELavg2  to 

%  get  the  tob  SEL  for  every  channel  of  that  Yokogawa.  It  assembles 

%  all  the  tob  SELs  for  a  particular  source  type  and  location  into  the 

%  output  matrix,  "results".  This  function  was  used  multiple  times  to  make 

%  tobSEL_all .mat 

%Requires:  Fcenter.mat,  MicLocations .mat 

%Calls:  SELavg2.m 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

srcType=  { ' C— 4 '  ' PC ' } ; 

srcLoc=  (1:4); 

scope= {' Trailer '  'Forest'}; 

cd  ('N:\Research\Acoustics\Camp  Ripley  Matlab  FormatX Spectrum  Analysis'); 
load  Fcenter.mat; 

cd  .  . 

load  MicLocations .mat; 

results { 1 , 1 } = ' SrcType ' ;  results { 1 , 2 } = ' SrcLoc ' ;  results { 1 , 3 } = 'MicLoc ' ; 
results {1,4}=' TotalSEL ' ; 
for  1=1:47 

results{l,i+4} =Fcenter ( 1 )  ; 

end 

count=2 ; 

for  i=l  %srctype 

for  j=l  %srcloc 

for  k=2  %scope 

if  strcmp (scope{ k} ,' Trailer ') ==1 

for  m=l : 1 6 

[Freq, total SEL, bands ] =SELavg2 (srcType (i) , srcLoc ( j ) , scope (k) ,m) ; 

for  p=l : length (totalSEL) 

results {count+p-1 , 1 } =srcType { i } ; 

results {count+p-1 , 2 } =srcLoc ( j  )  ; 

results { count+p-1 , 3 } =TrailerMicLocations {m} ; 

results {count+p-1 , 4 } =totalSEL (p) ; 

for  n=l : 47 

results {count+p-1 , n+4 } =bands  (p, n) ; 

end 

end 

results{count,  3} 
count=count+length (totalSEL)  ; 

end 


el seif  strcmp ( scope { k} , ' Forest ') ==1 
for  m=l:8 

[Freq, totalSEL, bands ] =SELavg2 (srcType (i) , srcLoc ( j ) , scope (k) ,m) ; 

for  p=l : length (totalSEL) 

results { count+p-1 , 1 } =srcType { i } ; 

results { count+p-1 , 2 } =srcLoc  ( j ) ; 

results { count+p-1 , 3 } =ForestMicLocations {m}; 

results {count+p-1 , 4 } =totalSEL (p) ; 

for  n=l : 47 

results {count+p-1 , n+4 } =bands (p, n) ; 

end 


end 
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results { count, 3 } 
count=count+length (totalSEL) ; 

end 


end 

end 

end 

end 
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function  [results] =SELperShot (srcstr, srcloc, config) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%INPUT :  srcstr  =  1=04,  2=PC 

%  srcloc  =  1  thru  4 

%  config  =  l=Parallel,  2=Perpendicular ,  3=Vertical,  4=all 

%OUTPUT:  "results"  matrix  containing  the  total  SEL  for  all  shots  matching 

%the  desired  input  values.  Columns  are  Mic  Location,  total  SEL,  and  shot  #, 

%respectively . 

%DESCRIPTION :  Filters  through  the  tobSEL_all  matrix  and  collects  all  the 

%shots  matching  the  desired  source  type,  source  location,  and  array  config. 

%These  are  plotted  to  show  the  total  SEL  for  each  mic  location  per  shot. 

%Requires :  tobSEL_all . mat 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

load  ('N:\Research\Acoustics\Camp  Ripley  Matlab  FormatX Spectrum  Analy- 
sis\tobSEL_all .mat ' ) ; 

if  srcstr==l 

srctype= ' C-4  '  ; 
elseif  srcstr==2 
srctype= 'PC'; 

end 

if  config==l 

array= ' Parallel ' ; 
elseif  config==2 

array= ' Perpendicular ' ; 
elseif  config==3 

array= ' Vertical ' ; 

end 

rescnt=0;  %rescnt  will  be  the  total  number  of  matches 

for  i=2 : length (tobSEL_all) 

if  (srcloc==tobSEL_all { i, 2 }  &&  strcmp ( srctype, tobSEL_all { i , 1 }) ==1  &&  ... 

(strcmp (array, tobSEL_all {1,3}) ==1 | | conf ig==4 ) ) 
rescnt=rescnt+l ; 

results { resent , 1 } =tobSEL_all { i , 5 } ;  %mic  loc  to  columnl 

results { resent , 2 } =tobSEL_all { i , 6 } ;  %totalSEL  to  column2 

results { resent , 3 } =tobSEL_all { i , 4 } ;  %shotnum  to  column3 

end 

end 

figure 
hold  on 
hand= [ ] ; 

for  k=l : length ( results ) 

if  (strcmp (results (k, 1 ),  'Ala ') ==1 I  I strcmp (results (k, 1 ),  'Alb ') ==1  I  I  ... 

strcmp (results (k, 1 ) ,  ' Ale ' ) ==1 I  I  strcmp ( results ( k, 1 ) ,  1  Aid1 ) ==1)  %Ala- 

Ald 

hand ( 1 ) =plot (results { k, 3 }, results { k, 2 },' bo ') ;  %hand  vector  will  be  used  for 

legend 

elseif  (strcmp (results (k, 1 ),  'A2a ') ==1 I  I strcmp (results (k, 1 ),  'A2b ') ==1  I  I  ... 

strcmp ( results ( k, 1 ) , ' A2c ' ) ==1 I  I  strcmp ( results ( k, 1 ) ,  ' A2 d ' ) == 1 )  %A2a-A2d 

hand (2 ) =plot (results { k, 3 }, results { k, 2 },' bs ') ; 

elseif  (strcmp (results (k, 1 ),  'A3a ') ==1 I  I strcmp (results (k, 1 ),  'A3b ') ==1  II  ••• 

strcmp (results (k, 1 ) , 1 A3c 1 ) ==1 [ I  strcmp ( results ( k, 1 ) , 1 A3d' ) ==1)  %A3a-A3d 

hand (3) =plot (results { k, 3 }, results { k, 2 },' bd ') ; 

elseif  (strcmp (results (k, 1 ),  'A4a ') ==1 I  I strcmp (results (k, 1 ),  'A4b ') ==1  ||  ... 

strcmp ( results ( k, 1 ) , 1 A4c 1 ) ==1 ] I  strcmp ( results ( k, 1 ) , 1 A4d' ) ==1 )  %A4a-A4d 
hand (4 ) =plot (results { k, 3 }, results { k, 2 },' bp ') ; 


%must  change  the  input  integers  to 
%match  the  values  in  the  cells  of  tobSEL_all 
%in  order  to  scan  it  later 


elseif  (strcmp (results (k, 1 ) , 'Ll') == 1 ) 

hand (5) =plot (results { k, 3 }, results { k, 2 } , 'mo ' ) ; 


%L1 
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elseif  (strcmp (results (k, 1 ),' L2 ') ==1 )  %L2 

hand ( 6) =plot (results { k, 3 }, results { k, 2 }, 'ms ') ; 

elseif  (strcmp (results (k, 1 ),' L3 ') ==1 )  %L3 

hand (7 ) =plot (results { k, 3 } , results { k, 2 } , 'md ' ) ; 

elseif  (strcmp (results (k, 1) ,  ' SI ') ==1)  %S1 

hand ( 8 ) =plot (results { k, 3 }, results { k,  2 },' go ')  ; 

elseif  (strcmp (results (k, 1) S2 ') ==1)  %S2 

hand ( 9) =plot (results { k, 3 }, results { k,  2 },' gs ')  ; 

elseif  (strcmp (results (k, 1) S3 ') ==1)  %S3 

hand ( 10) =plot (results { k, 3 } , results { k,  2 } ,  ' gd '  )  ; 

elseif  (strcmp (results (k, 1) S4 ') ==1)  %S4 

hand ( 11 ) =plot (results { k, 3 } , results { k, 2 } , ' gp ' ) ; 

elseif  (strcmp (results (k, 1) S5 ') ==1)  %S5 

hand ( 12 ) =plot (results { k, 3 }, results { k, 2 } , ' g+ ' ) ; 

end 

end 


hold  off 

xlabel ( ' Shot  Number') 
ylabel (' Total  SEL  (dB) ' ) 
set (gca, ' XTick ' , 1 : resent) 

title ( [sretype  '  -  BP'  num2str (sreloe)  '  :  '  array  '  Config']) 

leg¬ 
end  (hand,  ' A1 ' ,  ' A2 ' ,  ' A3 ' ,  ' A4 ' ,  'Ll',  ' L2 ’ ,  ' L3 ’ ,  ' SI ' ,  ’ S2 ’ ,  ’ S3 ' ,  ’ S4 ' ,  ’ S5 ' ,  1  Location ' ,  'NorthEas 
tOutside ' ) 
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function  [results] =  PeaksPerShot (srcin, srcloc, conf ig) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%INPUT :  srcstr  =  l=C-4,  2=PC 

%  srcloc  =  1  thru  4 

%  config  =  l=Parallel,  2=Perpendicular ,  3=Vertical,  4=all 

%OUTPUT:  "results"  matrix  containing  the  peak  level  for  all  shots  matching 

%the  desired  input  values.  Columns  are  Mic  Location,  peak  SPL,  and  shot  #, 
%respectively . 

%DESCRIPTION :  Filters  through  the  peaksTable  matrix  and  collects  all  the 

%shots  matching  the  desired  source  type,  source  location,  and  array  config. 
%These  are  plotted  to  show  the  peak  level  for  each  mic  location  per  shot. 

%Requires :  peaksTable . mat 

%Written  for  experiments  at  Camp  Ripley,  June  2006 
%by  Tim  Onder  (Timothy.E.Onder@erdc.usace.army.mil) 

Sfek  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

load  ('N:\Research\Acoustics\Camp  Ripley  Matlab  Format\Pk  vs  Dist  Analy- 
sis\peaksTable .mat ' ) ; 

rescnt=0; 

for  i=l : length (peaksTable ) 

if (peaksTable (i, 1) ==srcin  &&  peaksTable (i, 2) ==srcloc  &&  (con- 
fig==peaksTable (i, 3) I |config==4)) 
rescnt=rescnt+l ; 

results (resent, 1 ) =peaksTable (i, 4 ); %mic  loc  to  results  coll 
results (resent, 2 ) =peaksTable (i, 7 ); %peak (dB)  to  results  col2 
results (resent, 3) =peaksTable (i, 8 ); %shotnum  to  results  col3 

end 

end 

figure 
hold  on 
hand= [ ] ; 

for  k=l : length ( results ) 

if  ( results ( k, 1 ) >=1  &&  results ( k, 1 ) <5 )  %Ala-Ald 

hand ( 1 ) =plot (results (k, 3) , results (k, 2 ) ,  ' bo  '  ) ; 

elseif  (results (k, 1) >=5  &&  results (k, 1) <9)  %A2a-A2d 

hand (2) =plot (results (k, 3) , results  (k, 2) ,  'bs  '  ) ; 

elseif  (results (k, 1) >=9  &&  results (k, 1) <13)  %A3a-A3d 

hand (3) =plot (results (k, 3) , results (k,  2 )  ,  ' bd ' ) ; 

elseif  ( results (k, 1 ) >=13  &&  results (k, 1) <17)  %A4a-A4d 

hand (4 ) =plot (results (k, 3) , results (k, 2 ) , ' bp ' ) ; 

elseif  results ( k, 1 ) ==17  %L1 

hand (5) =plot (results (k, 3) , results  (k, 2 ) ,  'mo ' ) ; 

elseif  results (k, 1 ) ==18  %L2 

hand ( 6) =plot (results (k, 3) , results (k, 2 ) , 'ms ' ) ; 

elseif  results (k, 1) ==19  %L3 

hand (7 ) =plot (results (k, 3) , results (k, 2 ) , 'md ' ) ; 

elseif  results (k, 1) ==20  %S1 

hand ( 8 ) =plot (results (k, 3) , results (k, 2 )  ,  ' go ' )  ; 

elseif  results ( k, 1 ) ==2 1  %S2 

hand ( 9) =plot (results (k, 3) , results (k, 2 ) , ' gs ' ) ; 

elseif  results (k, 1) ==22  %S3 

hand ( 10) =plot (results (k, 3) , results (k, 2 )  ,  ' gd' ) ; 

elseif  results (k, 1) ==23  %S4 

hand ( 11 ) =plot (results (k, 3) , results (k,  2 )  ,  ' gp ' )  ; 


elseif  results (k, 1 ) ==24  %S5 

hand ( 12 ) =plot (results (k, 3) , results (k, 2 ) , ' g+ ' ) ; 
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end 

end 


hold  off 

xlabel ( ' Shot  Number') 
ylabel ( ' Peak  SPL  (dB) ' ) 
if  srcin==l 

source= ' C-4 ' ; 
elseif  srcin==2 
source= 'PC'; 

end 

if  conf ig==l 

arrayConf ig= ' Parallel '  ; 
elseif  conf ig==2 

arrayConf ig= ' Perpendicular '  ; 
elseif  conf ig==3 

arrayConfig= 'Vertical ' ; 
elseif  conf ig==4 

arrayConf ig= ' All ' ; 

end 

set (gca, ' XTick ' , 1 : resent) 

title ( [source  '  -  BP'  num2str (sreloe)  '  :  '  arrayConf ig  '  Config']) 
leg¬ 
end  (hand,  ' Al ' , ' A2 ' ,  ' A3 ' ,  ' A4 ' ,  'Ll',  ' L2 ' ,  ' L3 ' ,  ' SI ' ,  ' S2  '  ,  ' S3 ' ,  ' S4 ' ,  ' S5  '  ,  'Location' ,  'NorthEas 
tOutside ' ) 
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